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Abstract. A Gelfand-Tsetlin scheme of depth is a triangular array with 
m integers at level m, m = 1, . . . , N , subject to certain interlacing constraints. 
We study the ensemble of uniformly random Gelfand-Tsetlin schemes with 
arbitrary fixed Ath row. We obtain an explicit double contour integral ex- 
pression for the determinantal correlation kernel of this ensemble (and also of 
its g-deformation). 

This provides new tools for asymptotic analysis of uniformly random lozenge 
tilings of polygons on the triangular lattice; or, equivalently, of random stepped 
surfaces. We work with a class of polygons which allows arbitrarily large num- 
ber of sides. We show that the local limit behavior of random tilings (as 
all dimensions of the polygon grow) is directed by ergodic translation invari- 
ant Gibbs measures. The slopes of these measures coincide with the ones of 
tangent planes to the corresponding limit shapes described by Kenyon and Ok- 
ounkov |KO| . We also prove that at the edge of the limit shape, the asymptotic 
behavior of random tilings is given by the Airy process. 

In particular, our results cover the most investigated case of random boxed 
plane partitions (when the polygon is a hexagon). 



1. Introduction 

We study uniformly random tilings (by lozenges of three types) of polygons 
drawn on the triangular lattice with sides parallel to one of three possible directions. 
Equivalently, one can speak either of random stepped surfaces — that is, continuous 
3-dimensional surfaces glued out of 1 x 1 x 1 boxes with sides parallel to coordinate 
lines; or dimer coverings (= perfect matchings) of the hexagonal bipartite graph 
which is dual to the triangular lattice (see Fig. fTl). This model has received a 




Figure 1. A tiling of a polygon by lozenges {— rhombi) of three 
types, and its dimer interpretation. 

significant attention |CKPj , jKOj , jK2j . See also the lecture notes jK3| for a detailed 
exposition of the subject and more references. 
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The most well-studied particular case is the model of uniformly random tilings of 
the hexagon having sides a, 6, c, a, b, c (objects also referred to as uniformly random 
boxed plane partitions). Let us briefly indicate results obtained for the hexagon 
model, and how we generalize them in the present paper. A detailed description of 
the model and our results can be found in Sj2j 

1.1. Limit shapes for tilings of general polygons. The law of large numbers 
under the global scaling (when the polygon is fixed and the mesh of the lattice is 
going to zero) is well-studied for general polygons. In this limit, the random stepped 
surface (which is the same as the height function of the corresponding random dimer 
covering) approaches a non-random limit shape which can be obtained as a unique 
solution to a suitable variational problem. See |CLP| , [CKPj , and also |DMB| , [D] . 
This solution was described in [KO by means of the complex Burgers equation. 
For a wide class of polygons, the limit shape is an algebraic surface. 

Typically, such a limit shape forms facets, i.e., areas where asymptotically only 
one kind of lozenges is present. There is a curve in the plane (called the frozen 
boundary) that separates the liquid phase (where asymptotically one can see a 
random mixture of lozenges of all types) from the facets. This curve is also algebraic 
and is tangent to all sides of the polygon (or lines containing the sides) . See CLP] 
Fig. 2], [BG[ Fig. 5], |KOi Fig. 1] for illustrations of random tilings with small mesh 
where the limit shape and the frozen boundary are clearly seen. 

1.2. Asymptotics in the bulk. In the case of the hexagon it was established in 
[BKMM J and [Q (see also [Kl], |KOS] and [jT], [J2], [JN]) that locally near every 
point of the limit shape, the asymptotics of the random surface are directed by 
(uniquely defined |Shj ) ergodic translation invariant Gibbs measure on tilings of 
the whole plane. The slope of this measure coincides with that of the tangent plane 
to the corresponding limit shape of |K0) at the chosen point. 

The limiting translation invariant Gibbs measures are believed to be universal; 
they are present in random tiling models of, e.g., |0R1| . jOR2j . |BK| . |BMRTj . 
More general measures on tilings of the hexagon also exhibit the same local asymp- 
totics |BGR| . Furthermore, in |K2j such behavior was established for rather general 
random stepped surfaces with no asymptotically frozen zones. 

In the present paper we obtain the conjectural bulk limit asymptotics for uni- 



formly random tilings of a wide class of polygons (described in ^ 2.1 ). In particular, 
our approach allows polygons with arbitrarily many sides, which means that the 
frozen boundary curve and the limit shape surface can have an arbitrarily large 
degree. The slopes of the limiting Gibbs measures agree with the limit shapes of 
IKOI. 



1.3. Asymptotics at the edge. In the present paper we show that the asymp- 
totics of random tilings at the edge of the limit shape are directed by another 
universal law — the Airy process (introduced in [PS] in connection with a polynu- 
clear growth model). Until now even for the case of the hexagon only a partial 
result in this direction was established ^BKMM] via subtle analytical results about 
asymptotic properties of discrete orthogonal polynomials. Namely, the Airy process 
asymptotics were proven to hold only in one direction, i.e., for the statical ensemble 
of random particle configurations on a line, see |BKMMi Fig. 7]. We manage to 
obtain the (space-time) Airy asymptotics for tilings of polygons in our class. 
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It should be possible to go further and show (in a way similar to |0R2j . [BKj . 
and |BMRT] ) that at cusps of the limit shape the asymptotics of random tilings 
are governed by the Pearcey process. 

1.4. Connection with orthogonal polynomials. It is worth noting that almost 
all results cited above in §1.2| and §1.3| are obtained using in one form or another 
the technique of orthogonal polynomial ensembles. Because random tilings of poly- 
gons other than the hexagon lack a direct connection with orthogonal polynomial 
ensembles, no similar results about more general polygons were known. 

1.5. Technique of the present paper. We overcome the lack of orthogonal poly- 
nomials for more general polygons by establishing a new double contour integral 
formula (Theorem [T]) for the determinantal correlation kernel of random tilings. 
This provides necessary tools for the desired asymptotic analysis. 

We believe that our technique could also be applied to a more detailed study of 
various other models of random tilings such as, e.g., the ones considered recently 
in [BD], pm] and [NY2j . 

Acknowledgements. The author would like to thank Alexei Borodin for fruit- 
ful discussions, and Vadim Gorin for helpful comments. The work was partially 
supported by the RFBR-CNRS grants 10-01-93114 and 11-01-93105. 

2. Model and results 

2.1. Lozenge tilings of polygons. For convenience, we perform a simple afRne 
transform of lozenges which were present in Fig. [l] 

After this transform, the polygons will be drawn on the standard square grid with 
all sides parallel to one of the coordinate axes or the vector (—1, 1). We denote the 
horizontal and the vertical coordinates by x and n, respectively. 
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Figure 2. A polygon drawn on the ruled paper. In this example 
k = 4, and the polygon has 3k = 12 sides. 

We will restrict ourselves to polygons of a special kind as shown on Fig. [2] 
Every polygon P we consider can be parametrized by two integers N = 1,2, . . ., 
and fc = 2, 3, . . . (the polygon has 3fc sides) and by the (proper) half-integers 

Ai < Bi < A2 < B2 < ... < Ak < Bk, A„B, e Z' :=Z+ i, 
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subject to the condition J2i=i{^i ~ ^i) — The bottom side of P lies on the 
horizontal axis n = 0, and all the k—1 top sides lie on one and the same line n = N, 
so N is the height of P . 

The main object of the present paper is the uniform measure Pp on the set of all 
lozenge tilings of the polygon P by lozenges of three types (the transformed ones in 
(2.1 )). An example of such a tiling is given in Fig. |3] (we trivially extend the tiling 
to the whole strip < n < with N small triangles added on top). 




Figure 3. A tiling of P. Particles correspond to lozenges of one type. 

For fc = 2, our measure Pp becomes the uniform measure on all tilings of the 
hexagon with sides of lengths {A2 - Bi), (B2 - ^2)^2, (-Bi - Ai), {A2 - Bi), {B2 - 
^2)1/2, (Bi — Ai). The total number of tilings of such a hexagon is given by the 
classical MacMahon product formula, e.g., see [StJ §7.21]. Existing results about 
Pp in this case are discussed in the Introduction. 

There is also a product formula for the total number of tilings of a polygon P 
for arbitrary k, i.e., for the partition function of our model: 



^P = TT (2.2) 



where 



{xj^ < . . . < xf } = {Ai + i < Ai + I < . . . < Si - f < Si - i < (2.3; 

< A2 + ^ < ... < B2 - I < ... < Ak + I < ... < Bk - ^} 
are the positions of particles on the top line n — N, see Fig. |3] One way to obtain 



(2.2) is to interpret Zp as a dimension of a certain irreducible representation of the 
unitary group U{N), and then use the classical Weyl dimension formula ^Wj. See 
S|3]for more detail. 

2.2. Particle configurations and tiie correlation kernel. Another way to look 



at the number of tilings Zp (2.2) is provided by the dimer interpretation of our 
model (Fig. [T]) . Namely, Zp is equal to the determinant of the Kasteleyn matrix 
for the honeycomb graph (dual to the original triangular lattice) located inside P 
jKas] . In particular, this implies that the uniform measure Pp on tilings of P is 
determinantal (e.g., see |K3l Cor. 3]). There are several (very similar) ways to 
express the determinantal property of Pp, we choose the most convenient for us. 
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We put a particle in the center of every lozenge of type (see Fig. [sj) . Note that 
the coordinates of all such particles are integers. Thus, one sees a configuration of 
particles X := {x™ : m = I, . . . , N; j ^ I, . . . ,m} G Z^f^+i)/^ ^jth precisely 



m 



particles at the rnth horizontal line, m — 0, . . . , N. Because we have a tiling of P, 
these particles must satisfy the interlacing constraints: 

X™ 1 < x;"-i < X™ (2.4) 

(for all j's and m's for which these inequalities can be written out). Clearly, lozenge 



tilings of P and such interlacing arrays X with fixed top row as in (2.3) are in a 



bijective correspondence. These arrays are also in a bijection with Gelfand-Tsctlin 



schemes ([3.5) 



In this way, the measure Pp on tilings becomes a probability measure on inter- 
lacing particle arrays X = {x™}. We denote it also by Pp. 

Definition 2.1. Let (xi,ni), . . . , {xs,ng) be pairwise distinct points, e Z, 1 < 
Hi < N — 1. The correlation functions of the measure Pp are defined as 

Ps{xi,ni; . . . ;Xs,ns) := Pp (there is a particle of the random configuration {x™} 

at position (xj, Uj) for every j = 1, . . . , s). 

The dimer interpretation mentioned above implies that the measure Pp on inter- 
lacing particle arrays is determinantal. That is, there exists a function K{x, n; y, m) 
called the correlation kernel, such that 

Ps{xi,ni; . . .;xs,ns) = det[K{xi,ni;xj,nj)]lj^-^ (2.5) 

for any s and any collection of positions (xi, ni), . . . , {xs, Us). About determinantal 
point processes in general see the surveys jSosh] , |HKPVj , [B] . 

The first main result of the present paper is an explicit formula for the correlation 
kernel K of Pp in terms of double contour integrals: 

Theorem 1 (Correlation kernel). Fori <ni<N,l<n2< N—1, andxi,X2 E Z, 
the correlation kernel of the point process Pp has the forn^ 

(xi - a;2 -f l)„i_„2-i {N-ni)l 



K{xi,ni;X2,n2) = -ln.2<nilx2<xi- 



(ni-n2-l)! (A^-n2-l)!' 



1 I I {z-X2 + l)w-„2-l 1 A + 5 ^ '^)i3.-A. 



/ — TtIat ,1 7/1— 



(27ri)2 / / {w-Xi)N-n,+i w - z ^J- (A, + ^ - z)b^-a,' 

{z} {w} 

(2.6) 

The contours in z and w are positively (counter-clockwise) oriented and do not 
intersect. The contour {z} encircles the integer points X2, X2 + 1, ■ ■ ■ , Bk — | in z 
and only them (i.e., does not contain a;2 — 1, 2:2 — 2, . . .). The contour {w} contains 
{2} and all the points a;i, Xi — 1, . . . , Xi — {N — ni). 

We prove Theorem [T] in §f|4]-[5] using an Eynard-Mehta type theorem which al- 
lows the number of particles to vary. Its statement can be found in, e.g., pBj §4]. 
However, the application of that theorem is not straightforward. Namely, to get 
the result about the uniform measure Pp, we first consider its g-deformation ^Pp 
obtained by weighting every tiling proportional to q"°\ where vol is the volume 



^Here and below !{...} denotes the indicator of a set, and (y)m '■= y{y -I- 1) • ■ • (s/ + ™ ^ 1), 
m = 1, 2, . . . (with {y)o ■= 1) is the Pochhammer symbol. 
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under the corresponding stepped surface. Such deformed measures were consid- 
ered in |K0| and (as a particular case) in |BGR| . The correlation kernel qK of 



is given by a double contour integral expression similar to (2.6), but with a 



q-hypergeometric function inside the integral (see Theorem 4.1 1. 

Then in a limit as 9 1 we obtain the formula for the kernel K of Theorem [T] 
The complexity of the expression for qK is the reason why restrict our asymptotic 
analysis to the case q = 1. The disappearance of the g-hypergeometric function in 
the g — > 1 limit may seem rather surprising. A possible explanation can be given in 



the course of computations (see Remark 5.8 and Proposition 5.9): the limit 1 



kills many additional negligible terms, and this significantly simplifies the post-limit 
{q = 1) formulas. 

Remark 2.2 (Kasteleyn matrix). Using Theorem [T] in Sj6]we manage to express 
the inverse Kasteleyn matrix for the honeycomb graph inside P (see Fig. [ij right) 



through the correlation kernel K (2.6). Similar results for other models can be 



found in |BF2l §5.1] and |BGR[ §7.2] (the latter model includes the hexagon case). 

Remarks 2.3 (Degeneration to some known kernels). 1. A similar kernel for the 
continuous Gelfand-Tsetlin patterns was obtained (using similar methods) recently 
in |Met[ Prop. 2.4]. The model of |Metj describes distributions of eigenvalue minor 
processes of random Hermitian N x N matrices. It can be shown that the kernel of 
[Met) is a limit of ours when one: 1) embeds every horizontal line n = 0, . . . ,n = N 
(see Fig. jsj) into M as Z 1— > ^Z; 2) scales the positions of particles as = 
with new continuous a:"'s; and then lets ^ — > 0. Note that in this continuous limit 
the connection with random tilings is lost. 

2. In |BK| . a correlation kernel for so-called ergodic measures on (infinite) 
Gelfand-Tsetlin schemes was obtained. Equivalently, one can speak about random 
tilings of the whole upper half plane when the tiling spreads all the way upwards 
(N = 00 on Fig. |3| . These measures are limits (in some sense) of the measures 
Pp as — > 00. This limit transition is related to classification of characters of 
the infinite-dimensional unitary group. See the recent paper |BO) and references 
therein for more detail on the subject. It can be shown that the kernel of [BKj is 
a limit of our kernel K in this regime. 

2.3. Asymptotics in the bulk. We consider A^ — > 00 asymptotics as all dimen- 
sions of the polygon P(A^) grow. We assume that the parameters Ai{N), Bi{N) of 
P(A^) are scaled as 

Ai{N) ^ [aiN] + const, B,{N) ^ [hN] + const {i = l,...,k). (2.7) 

Here ai < bi < ... < < bk are new continuous parameters which satisfy 
Si=i(^i — tti) — 1. The constants above are bounded uniformly in A^ and are 
needed to ensure that Ai{N),B,{N) e Z' and J2k=liB^{N) - Ai{N)) = A^. 

One can then rescale the growing polygons P(A) by A^~^ in both directions, 
so they will approach some fixed polygon V drawn on the new (%, rj) coordinate 
plane. The polygon V is parametrized by {ai,6i}fL]^ in the same way as it was for 
P(A^) and {Ai{N), B,{N)}'^^j^, and V is located inside the strip < r] < 1. See the 
polygon on Fig. |4] (we preserve the proportions of the polygon on Fig. [2j|3]). 

Alternatively, one could say that we consider tilings of a fixed polygon V with 
refining mesh. In terms of particle configurations, this asymptotic regime means 
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that the N particles on the top row form k macroscopic clusters in which the 
particles are densely packed. 



fl4 64 




Figure 4. The limiting polygon V on the {x,v) plane and the 
frozen boundary curve. The segments [ai,bi] indicate densely 
packed clusters of particles on the top row. 



Let us discuss the limiting object which will describe local asymptotics of random 
tilings. 



Definition 2.4 (Incomplete beta kernel [ORT]). LetneC\ {0, 1}, 3f7 > 0, be a 
parameter called the complex slope. The incomplete beta kernel Bf2(m,^) is defined 
as 



Bn{m,l) := 



1 



(1 - u)"'u''-'du, 



m, / e Z, 



27ri JO 

where the path of integration crosses (0, 1) for m > and (— oo, 0) for to < 0. 



This kernel was introduced in [ORl] in connection with local asymptotics of 
random plane partitions. As was shown in |Shj and [KOSj . for every complex 
slope there is a unique ergodic translation invariant Gibbs measure on tilings of 
the whole plane, and this is the determinantal point process corresponding to the 
kernel Bf2(TOi — TO2, ^1 — /2) (as in (2.5), the correlation kernel depends on two points 



(toi, li), (to2, ^2) in the plane). The complex slope 51 defines proportions of lozenges 
of various types Ai,Pn,p^, P^+Pn+P-^ = 1 (seen in a large box under this measure) 
as follows. For Qfl > 0, p<^,Pa,p-z^ are proportional to angles of the triangle in the 




Figure 5. Lozenge densities defined by a complex slope fl, 317 > 0. 
complex plane with vertices 0, fl, 1 (see Fig. [5]). For G M, 7^ 0, 1, the Gibbs 
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measure is supported by a single configuration with all lozenges of one type. This 
may be seen as a degeneration of Fig. [5] when fl approaches the real line. 

We note that the cases 17 = 0, 1, or oo in our model correspond to points where 



two different frozen facets meet (see Proposition 2.7). At these points there is no 
translation invariant local asymptotics. 

Theorem 2 (Asymptotics in the bulk). Let (XtV) encode a global position in the 
scaling limit (see Fig. [^. As N — ^ oo, locally around (Xi*?); the measure Pp(Ar) 
on random tilings converges to an ergodic translation invariant Gibbs measure on 
tilings of the ivhole plane with some complex slope = r2(x, rf). In terms of correla- 
tion functions, this means that for any points {xi{N),ni{N)), . . . , {xsiN),ns{N)) 
such that 

Xi{N)/N-^X, Ui/N rj, N oo {i = l,...,s) 
while the differences Xi — Xj , — rij G Z stabilize, we have the convergence 
p,{xi{N),ni(N); . . . ; x,(iV), n,(iV)) det[Bo(n, - n^ixj - a;,)],"^^!. 
This theorem is proved in ^ using the double contour integral representation 



(2.6) for the correlation kernel K. The proof uses the saddle point analysis and 



mainly foUows the lines of [Ok], [ORl] . [0R2] . and [BK] 

2.4. Complex Burgers equation, limit shape, and frozen boundary. Let 

us describe how the complex slope D,{x,ri) at a global point {x,v) (Theorem [2| is 
characterized. 

There is an (open) domain I? inside the polygon V (the so-called liquid region) 
such that for every {x,il) G 2^ we have Qn,{x,ri) > 0. This means that inside T) all 
types of lozenges are asymptotically present (in proportions determined by r2(X; 'y); 



see S2.3) 



On the other hand, everywhere in 7^ \ I?, the random configuration is asymp- 
totically frozen, i.e., it contains only one type of lozenges. The liquid region is 
separated from frozen ones by a curve dT> called the frozen boundary (see Fig. |4]). 
For our model, dV is an algebraic curve of degree k inscribed in 7^: it is tangent 
to all sides of the polygon V (or lines containing them) . The curve dV has k — 2 
turning points (corresponding to cusps of the limit shape) . In the k — 2 case when 
■p is a hexagon, dT> reduces to the ellipse inscribed in V. 

Proposition 2.5 (The complex slope n{x,ri)). For {x,v) ^ slope fl{x,ri) 

is the only solution in the upper half plane of the following algebraic equation^ 

k k 

^ ■ n -X + l~ri)^- (a,: " x)) = H ((^' - X + - v)^ - {h - x))- (2-8) 

1=1 1=1 

The slope Q.{x,ri) also satisfies a differential equation (called the complex Burgers 
equation |KOj ); 

mx,.)^ = -(l-^^(x,.))^. (2.9) 



^This equation has degree k + 1, but it always has the root f2 = 1, so in fact Q satisfies a 
degree k equation. There are other (s ometi mes more suitab le) c omplex par ameters of the limit 



shape and the frozen boundary — Wc i 2.10 i and T (Remark 2.8 I. See also §[ 7.6 7.7 
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In particular, this implies that the complex slope Q{x, v) of the limiting ergodic 
translation invariant Gibbs measure describing the local asymptotics around the 
global position (XtV) (Theorem [2]) coincides with the complex slope of the tangent 
plane to the limit shape of [CKP,, [KOj at the point (x,?/). This means that the 
normal vector to the limit shape at every point (x, rj) Cz V\ dT> looks as {p^, Pm P'^) 
(in the coordinates shown on Fig. [6]) . On the frozen boundary dV the limit shape 
surface does not have a normal vector. This interpretation is the reason why we 
call the parameter the complex slope. 




Figure 6. 3-dimensional coordinate system for a random stepped surface. 

Our analysis also produces an explicit rational parametrization of the frozen 
boundary dT> (which, in particular, is used to draw the curve on Fig. |4j see also 
Fig. 13 1. For {x,v) G 9T>, the two complex roots 51,^2 of (2.8) coincide. Thus, one 
could take 57 e M as a parameter for the curve. However, it is more convenient to 
use another parameter Wc, where 

Wc:=X+— i 7^, ^= — , • (2.10 

Note that is monotonically increasing in w^. 

Proposition 2.6 (Frozen boundary). The frozen boundary curve can be paramet- 
rized as 

I ^ n(wc)-i n(w,)(s(wc)-n(w,) + 2)-i 

X(wc) = + ; ry(w,)^ , (2.11) 

L(Wc) li(Wc)Z.(Wc) 

where 

Ui^c) :^ n' , - , - (2-12) 

-«-J-i=l Wc — Oi ^ — '1=1 \Wc~ bi \Nc — Oi/ 

with parameter — oo < w^, < oo. As increases, the frozen boundary is passed in 
the clockwise direction (so that V is to the right of dT>). 

We can also identify the tangent points on the frozen boundary: 

Proposition 2.7 (Tangent points). The tangent vector (x(wc), ?7(wc)) to the frozen 
boundary curve has slope 

X(Wc) Wc - x(wc) . . , . 

—, — V = ~i f — r =■ ' (.*=)• (2-13) 

There are sides of V of three directions to which the curve dT> is tangent: 
(S) The values w^, = a^, i — I, . . . ,k, correspond to points of dT> where it is tangent 
to sides of the polygon V parallel to the vector (—1, 1). Here ^{oi) = oo. 
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(V) When Wc ~ bi, i = 1, . . . ,k, the frozen boundary dV is tangent to vertical sides 
ofV. We have fl{b,) = 0. 

(H) At points \Nc — hi, . . . , hk, where il{hi) = 1, the frozen boundary dT> is tangent 
to horizontal sides of V . This includes the case hi — —oo where the curve is 
tangent to the bottom horizontal side ofV . Clearly, hi £ [bi-i, Oi), i — 2, . . . ,k. 

Finally, in each segment (/i2, h^), . . . , {hk-i, hk) there is one turning point of the 

frozen boundary (corresponding to a cusp of the limit shape). 

See also Fig. 15 for graphs of T{vJc) and fl{wc) along the frozen boundary curve. 
Rema rk 2. 8. Let us extend the definition of T (2.13) inside the liquid region 



(using (2.10)) as T{x,v) 



Wc(x.'?)-X 



Observe that the complex Burgers equation 



(2.9) can be rewritten in the following form inside V: 

dx dr/ 



T{x,v)- 



(2.14) 



The (intuitive) difference in sign between (2.13) and (2.14) indicates that the equa- 
tion (2.14) (and thus also ( |2.9| )) which holds inside the liquid region V does not 
survive the limit transition to the boundary of V. 

2.5. Asymptotics at the edge of the limit shape. Here we explain our results 
about the asymptotics of random tilings when the global position (%, 77) we are 
looking at is located on the frozen boundary. See S|8]for precise formulations. 

An intuitive way of thinking about the edge asymptotics governed by the Airy 
process (introduced in [PSJ) is to consider nonintersecting path ensembles. Namely, 
the frozen boundary at the chosen point (x, rj) can have one of three directions 
(named according to Proposition 2.7), see Fig. [7j Note that this is the same as to 
classify frozen regions by the type of lozenges they are built of. Close to a point at 




Figure 7. VH (left), HS (center), and SV (right) directions of the 
frozen boundary, and the corresponding nonintersecting path en- 
sembles. 

every direction of the frozen boundary, let us look at the corresponding ensemble 
of nonintersecting paths constructed using our random tiling (according to three 
cases on Fig. [t]) . We consider distribution of several such paths starting with the 
outer ones. That is, we look at first appearances of lozenges of two new types 
when the frozen region turns into a liquid one. These first nonintersecting paths 
concentrate around the corresponding part of the frozen boundary (when iV — > 00, 
after a global scaling by N^'^ in both directions). 
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The Airy-type asymptotics we establish for these paths means that at every 
given position (x, f?) G dV (we assume that (x, 77) is neither a tangent nor a turning 
point), the fluctuations of paths in direction tangent to dV are of order A^^/^, and 
in normal direction they have order A^^/^ (without a scaling). This can be restated 



in terms of our point processes Pp(Ar) as follows (see Theorem 8.1 for a full and 
detailed statement): 

Theorem 3 (Asymptotics at the edge). After a proper rescaling (and additional 
replacement of particles by holes and vice versa for the SV part of dV), the corre- 
lation functions of the point process Pp(Ar) converge to those of the Airy process. 
The latter are given by minors of the extended Airy kernel: 

, f/o" e-"(-i--^)Ai(ai + u)At{a2 + u)du, if > T2; 

(2.15) 

Here Ai{x) is the Airy function 

Ai{x) = — / e'*'/3+»xt^^_ (2.16) 
27r J -00 

We establish this Theorem also using saddle point analysis in a way similar to 
[0R2j . and BK . 

The rest of the paper is devoted to proving Theorems [T] [2] and [3j 

3. Combinatorics of the model and random Gelfand-Tsetlin schemes 
We will argue in terms of measures on interlacing integer (particle) arrays X — 



{x™ : m = 1, . . . , N; j = 1, . . . , m} (S 2.2 ) with an arbitrary fixed top row. In this 
section we explain the connection of such arrays with Gelfand-Tsetlin schemes and 



discuss the partition function of our model (2.2) and of its (/-deformation 



3.1. Signatures and branching of Laurent-Schur polynomials. By a signa- 
ture of length A^ we will mean a nonincreasing A^-tuple of integers A = (Ai > . . . > 
\n) G Let GTn denote the set of all signatures of length A^. Because GT^v 
parametrizes irreducible representations of the unitary group U{N) jWj, in the lit- 
erature signatures are also referred to as highest weights. The irreducible characters 
of U{N) are given by the Laurent-Schur polynomials jWj . [M] 



det 



u, 



sx{ui, . . . ,Ujv) = / ' N-,.M^^ ^ ' AgGTtv- (3.1) 



det[ 



Hj = l 



Here Ui, . . . , uj^ are eigenvalues of a unitary matrix U G U{N). Every s\{ui, . . . ,upf) 
is a rational function which is symmetric in ui, . . . ^uj^. More precisely, it is a 
homogeneous Laurent polynomial (i.e., a polynomial in uf^^...,u^^) of degree 
|A| := Ai -I- . . . -I- Ajv G Z (this number is not necessary nonnegative). Note that the 
denominator in (|3.1| is simply the Vandermonde determinant 



V{u,, . . . , ^i^.) det[u^-%^, =11,^ ^ ^^i)- (3-2) 

The branching of the Laurent-Schur polynomials reflects the branching of the 
irreducible characters of unitary groups: 

Sa(ui, . . . ,MAr_i;UAr = 1) = S^(mi, . . . , UjV-l), AgGTat, 
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where /U ^ A means interlacing of signatures: 

Al > /^l > A2 > . . . > Aat-i > > Aat. 

In fact, more can be said: 

s\{ui, . . . ,un^i;un) s^(ui, . . . ,MAr_i)u'^'~'^'. (3.3) 

^ — 'fj,: fj.-<\ 



Continuing expansion (3.3) for itjv-i, 'Ujv-2j ■ ■ ■ we arrive at 



Definition 3.1. A Gelfand-Tsetlin scheme of depth N is an interlacing sequence 
of signatures 

0^A(i)^A(2)...^A(^-i)^AW. 

One can also view this as a triangular array of integers {A^'"''} G 2^*^^+^^/^ with 
the following interlacing constraints: 

UN) j.(N) 



polynomials 



From (3.3) we get the following combinatorial formula for the Laurent-Schur 



M Ch. I, (5.12)]: 



sx[ui, . . . ,un) — 2^ U2 •■•Ujv ■ (.'J-'^i 

0^A(i)^...^A(")=A 

Here the sum is taken over all Gelfand-Tsetlin schemes with fixed top row A'^-* = A. 
We will use this identity in the present section to compute the partition function 
in our model and in its q-deformation. 

3.2. Volume of a Gelfand-Tsetlin scheme. In any signature A = (Ai, . . . , Aat), 
one can separate the positive and the negative components, and write it as a pair 
of partitions {— Young diagrams M, Ch. I.I]) A = (A+, A^): 

A = {X'l , . . . , A^ , 0, . . . , 0, — A^_ , . . . , — A]^ ), 

where A± > . . . > A±± > 0. Clearly, |A| |A+| - |A-|. 

Each Gelfand-Tsetlin scheme of depth TV with A^th row can 
be viewed as a pair of sequences of Young diagrams (A(™))±, m = l,...,N, where 
each (A^"*-*)^ is obtained from (A^™^^-*)* by adding a horizontal strip [Ml Ch. I.l]. 
We can make this into a pair of 3-dimensional Young diagrams (plane partitions) 
by placing N — m boxes of dimensions 1 x 1 x 1 on top of each 2-dimensional 
horizontal strip (A(™))±/(A(™-i))±. Let vol=^ = vol=^(0 -< X^^'> ^ ... ^ A^^)) be 
volumes of the resulting two 3-dimensional Young diagrams. 

More formally, 

vol± := . - ™)(I(A^™^)^| - |(A("-i))±|) = Y.""'] |(A("))±| 
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(by agreement, KA^^-*)*! = |0| = 0). Let us call the quantity 

vol :^vol+-vor =y^"V"^l (iV-m)(|A(")|- |A("-i)|) (3.5) 

the volume of the Gelfand-Tsetlin scheme -< A^^' ^ ... ^ xi^-i) ^ ^(w) ^ x. 

3.3. Measure g^™' on Gelfand-Tsetlin schemes. Everywhere in the paper we 
assume that < q < 1. 

Let us fix any signature v E GT ^ . Wc consider the following family of probability 
measures (depending on the parameter q) on the set of all Gelfand-Tsetlin schemes 
of depth N with fixed A^th row i^: 

-vol(0^i/<i)^...^iy("-i'^i/) 



-< V 



(N-l) 



(3.6) 



-vol 



We will sometimes abbreviate qVN,v as q' 

From the combinatorial formula (3.4), we readily conclude that the partition 
function is 

ry _ I 1-N -1 



,q -,1). (3.7) 
This specialization of Schur polynomials is known [M, Ch. I, §3, Ex. 1]: 



'9 '^^ <? 



(3.8) 



where V is the Vandermonde determinant (3.2). 



3.4. Uniform measure. Taking the g t 1 limit of qf'N,u, we arrive at the uniform 
measure on all Gelfand-Tsetlin schemes of depth N with the fixed A^th row v e 

PAr,,(0 < 1/(1) ^ . . . ^ -< v) := (3.9) 



It is known that for Vn,v the partition function has the form (cf. (3.7)-(3.8)) 



s.(l,...^) 

N 



T^(-l,. 



. , i^AT - A^) 



(3.10) 



From the branching rule (3.3) one can in fact conclude that Z^^v is equal to the 



dimension of an irreducible representation of U{N) corresponding to the signature 
pGGTn [Wj. 

3.5. Connection to random tilings. Gelfand-Tsetlin schemes -< z/'^^-' ^ . . . ^ 
j^(N-i) ^ ^(N) bijectively correspond to interlacing integer arrays X — {x™} (dis- 



cussed in \ 2.2 in connection with tilings) by means of a simple transformation: 



TO=l,...,iV, j = 1, 



The interlacing conditions for Gelfand-Tsetlin schemes of Definition |3.1| are clearly 



translated into the constraints (2.4) which are satisfied by X. 
We will think of the measures , 



and fN,u introduced in §^ 3.3 3.4 



as measures 



on such interlacing arrays (and do not change the notation for the measures). 



In the language of ^2.1 the uniform measure Pat^^ clearly becomes the corre- 



sponding uniform measure on all tilings of the polygon obtained by fixing the top 
row particles (see Fig. |3|) at positions xj^ = Vj — j, j = 1, ... ,7V. In particular. 



formula (2.2) for the partition function of random tilings is the same as (3.10) 
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Now let us consider g-deformed measures. For random tilings, one way to define 
volume vol under the corresponding stepped surface is to postulate that the tiling 
with zero volume is given on Fig. [8j Then the surface corresponding to every other 
tiling is obtained from the one on Fig. [8] by adding vol boxes of dimensions 1x1x1. 
It is readily checked that thus defined volume vol under the stepped surface is (up to 
a constant summand) equal to (—vol), where vol is the volume of the corresponding 

This implies that weighting of tilings proporti ona l 
is the same as considering the measure qfj^i, (3.6) 



Gelfand-Tsetlin scheme (^3.2 
to g""' (as mentioned in [2.2 
on interlacing arrays. 



—vol 



4. Correlation kernel for the measure q 

4.1. Formulation of the result. In this section we establish a q-deformation of 
Theorem [Tj that is, compute the correlation kernel qK of the measure qVn,!^ on 
interlacing arrays with fixed A^th row ([j3]). 

To formulate the result, we will need the standard notation of the g-Pochhammer 
symbol 

ia;q)k:^Yl'''a-aq') (fc = l,2,...), (a;(z)o:=l, 
and the g-hypergeometric function 

2<t)i{a,b;c\ q;z) := 



{c;q)t iq;q)i 



Theorem 4.1. The correlation kernel qK of the measure qPN.i^ on interlacing 
particle arrays {x™ : m = 1, . . . , N, j — 1, . . . , ni} with fixed top row xj^ = i^j — j 
(j = 1, . . . , N ) is given for I < ni < N , 1 < n2 < N ~ 1, and xi, X2 € Z by 

qK{x,,n,-X2,n2) = -ln,<n,l.,<.,g"^("^-"^) ^'^'^ W' (4-1) 

\q\ q)ni-n2-l 

{q^^^;q^^)N-m / / ^"^(^i^^^)^"^ 



(27ri)2 7 J w 

C,(a;2— Ki) c(oo) 



X ,0i(g-\g--i;g^-i I q-';w-') '''' . . . ^ \{ 



^.izq^-^-+^^;q)M-n,-l ^ w - q'^'-^-^^ 



iq;q)N-n2-i z-q^- ' 

The contours in z and w are counter-clockwise and do not intersect. The con- 
tour (tq{x2 — xi) in z encircles the points g^^-Ki ^ gX2-xi+i^ ^ _ ^ ^^^\-\-xx ^ ^^^^ 
only them (i.e., does not contain g^2-a;i-i^ ^2:2-2:1-2^ ^ ^ rp^^ contour c(oo) in w 
contains ^q{x2 — a^i) and has a sufficiently large radius. 

We prove Theorem |4.1| in the rest of this section. 
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4.2. Writing the measure q ™' as a product of determinants. In the proof of 
Theorem |4.1| we follow an Eynard-Mehta type formalism which allows the number 
of particles to vary. See iB] Thm. 4.4] for a detailed formulation of the statement 
we use (we will recall main formulas below in { 4 



be found in, e.g. [BFl] . [FM §4.4], see also [B 



4.3 U.5 and §4.6^. Its proofs can 



'PS[ Lemma 3.4]. The application 



of this formalism to our situation is in some aspects similar to |BKI . 

The first step in applying the formalism is to write the measure ^Pat.^ on in- 
terlacing arrays {x™} as a product of determinants. Since we allow the number of 
particles to vary, we add to every mth row of the array x™ < . . . < x™ a virtual 
particle "X-m+i — virt. Informally, one can think that virt = —oo. 

Let us denote for m = 1, . . . , A^: 

^,„(x, y) g(™-i)(^-"h,<, + g^^-^^n^*, x G {vtrt} U Z, y € Z. (4.2) 
The following lemma is well-known (e.g., see [BK, §3]): 

Lemma 4.2. For any array of integers {x™ : m — 1, . . . , N; j = 1, . . . , to}, we 
have (with the agreement that x™_)_]^ = virt) 

1 v'Ml". - l'<""'L(JV-i)|x"| 



if the array is interlacing (i.e., satisfies 1^2.4^ ), and otherwise. Here and below 
we denote jx'"] := x^" . . . x™. 

We need to define one more object which we will use: 

Definition 4.3. Let q\/{v) denote the x Vandermonde matrix [{q'^"^^^)^ j^i. 
Clearly, det[gV(z^)] = V{q'^^^^ , . . . , q^'^~^), where V is the Vandermonde determi- 
nant (3.2 1. 

Let gV(i/)"i be the inverse of that Vandermonde matrix. Define the following 
functions in x G Z: 

Aix 1 N) := ^^.^J,V(j.)-i],,l,.,^_,. (4.3) 



N.i 



Proposition 4.4 (Measure q as a product of determinants). The measure ql 
on arrays of integers X = {x™} G (mith x™_|_]^ — virt) can be written in 

the following form: 

,PjV,.(X) = qhmN+l)i2N+l)y(^^-l^ q^N^^ (4 4) 

N 

X det[V.(xf 1 7V)]5^i n det[^^{xr\x])]l,=^. 

n=l 

The proposition in particular claims that if the array is not interlacing or its 
top row x;j^, . . . , does not coincide with (i^i — 1, ... , ly^ — N), then the measure 
assigned to such array by the right-hand side of (4.4) is zero. 



Proof. By Lemma |4.2[ the product of determinants over n — 1, . . . , A'^ in the right- 
hand side of (4.4 1 imposes the interlacing constraints on {x™}. It also provides the 
desired weighting proportional to because jx™] — ]i/(™)] — ('"^^) (see ^3.5|. 

Then an obvious (but crucial; see Remark 4.9 below) observation is that for any 
integers yi > . . . > y^ one has 
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This allows to take the factor 



from the partition function qZ^^^ 



y(<j''l-l,...,g''iV-«) 

of qlPjv.i/ (see ^3.3) and make it into the determinant det[i/'i(xj^ | N)]fj^-f^. This 
imposes the desired condition that the top row of the array {x™} is fixed: xj^ = 
i^j -j,j = '^,---,N. 

The factor V{q~^, also comes from the partition function qZjsi^u- Fi- 
nally, it is readily checked that the power of q in front in (4.4) is correct. This 
concludes the proof. □ 



4.3. Convolutions and the "Gram matrix" . In § §4.3 4.5 using the expression 
of the measure qPjv.i/ as a product of determinants (Proposition 4.4 1, we compute 
quantities which are needed for the Eynard-Mehta type formalism of [51 Thm. 4.4] . 

Formula (4.9) for the kernel involves certain convolutions of the functions ipn 
(4.2). The computation of these convolutions is done similarly to [BK. §3]. Namely, 
define for ni, 712 — 1, . . . , N and a;, y € Z: 

The "*" notation means convolution, i.e., summation over an intermediate variable. 
See [B] §4] for a precise definition. 

Lemma 4.5. We have for x,y E Z; 

(9; 9)Ti2-ni-l 

Proof. Denote Fn{z) :— ■ For any x, y € Z, we clearly have 

^"("'^) = iijj,|^, 

Then the convolutions take the form 

1 



(ni,n2) 



dz „ , , 



(ni,ri2) 



dz 



= 1„,<„,V-('?"S-- (4-5) 
where (to G Z) are the complete homogeneous symmetric polynomials, /iq = 1, 
h^i = /i_2 = . . . = [Ml Ch. 1.2]. We have used their generating function 

00 L -. 

Y.h^{y„...,yL)t"^^l[^—-^. (4.6) 

m=0 r=l 

Since the h^s are particular cases of the Schur polynomials, = lm>oS(m) [Ml 
Ch. 1.3], we can write using (3.8) for every L > 1: 



^m(i,g, 



) = 1 



r?j>0" 



V{q"^-\q-^, 



Viq-\...,q-L) 



^m>0 



L-1 



q 



,m+l 



(.1 



q)L-i 



{q;q)L-i 

Since is a homogeneous symmetric polynomial of degree m, we see that the 



above expression for hm{^,q. 



,q 



L-l 



) together with (4.5) implies the claim. □ 



ASYMPTOTICS OF RANDOM LOZENGE TILINGS VIA GELFAND-TSETLIN SCHEMES 17 



We also need the following convolutions with one of the variables being virtual 
(here ni, = 1, . . . , and y G Z): 

Lemma 4.6. We have for y G Z; 

112— ni 
r=l 

(of course, for ni — n2 the empty product is interpreted as 1). 

Proof. For ni = n2 the claim is trivial. For rii > n2 from the proof of Lemma |4.5| 
we have: 



— '-ni<n2H / , 'i'y—x\<d, ■ ■ ■ i H J 

X 



xez 

(rii-l)i/ 



lni<n2 



Li the last summation we used (4.6) with t = 1. Note that in this case the series 
converges because < g < 1. □ 



Remark 4.7. The need for a converging series in the proof of Lemma 4.6 is the 
reason why one cannot directly use the Eynard-Mehta type formalism to compute 
the correlation kernel K of the uniform measure fN,v Thus, to get a formula for 
g = 1, we first need to deal with the case < g < 1, and then take the limit <? t 1- 
In fact, a similar problem occurs in computations in |BK1 §3]. 

The "Gram matrix" for our model is defined as 

G,ni := {iprn * ...(pN * V'K' I N)){virt), m,l = l,...,N. 

Lemma 4.8. We have 

N—m 

Gml — h=N-m+l ■ W_ ^ l"^) ^ ■ 



Proof. Let us write (by Lemma 4.6 1 

N-m 

{^„,*...^N*M-\N)){virt)= W (l-g'-)-'I]9^""'^V/(y|^)- 
Using the definition of ipi ( |4.3[ ), we can readily simplify the sum over y above: 

EN 
. Aq'^i'^) ]l,3 ■ [9V(l/)]j-7V-m+l 
J = l 

= l/=Af-m+l. 

This concludes the proof. □ 
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Remark 4.9. We see that the "Gram matrix" G — [Grni]m i=i turns out to be 
diagonaL This means that we readily know its inverse which enters the Eynard- 
Mehta type formula for the kernel qK ( 4.9 ). In fact, in various determinantal models 
inverting the matrix G is the main technical difficulty in obtaining an explicit 



formula for the correlation kernel. For our measure 



functions ^l)i{x \ N) in Proposition 
matrix" . 



4.4 



it is the use of the 



that allows to obtain a diagonal "Gram 



4.4. Inverse Vandermonde matrix and double contour integrals. Here we 
show how the elements of the inverse Vandermonde matrix can be written as double 
contour integrals. In fact, this is the reason of appearance of the double contour 
integral expression for qK in (4.1). 



Proposition 4.10. For every i,j — 1, 



1 

(27ri)2 



dz 



w 



N , we have 

dw 1 

w — z 



N+l-i , 



N 

n 



w — q 
z - a" 



(4.7) 



c(oo) 



The contour c{q'^^^^) in z is counter-clockwise, sufficiently small, and goes around 
q'^^^^ (it does not include any other poles in z). The counter-clockwise contour 
c(oo) in w contains c{q'^^^'->) (without intersecting it) and has a sufficiently large 
radius. 

Proof. Using the elementary symmetric polynomials Cm (e.g., see [Ml Ch. 1.2] for 
definition), one can write for i, j — 1, . . . , N: 



,V(^) 



■ ,9 



(4.8) 



(hat means the absence of q'^^^^). Indeed, this formula follows from the fact that 
every cofactor of the Vandermonde matrix gV(j/) can be expressed through the 
numerator in the right-hand side of (3.1) with A of the form (1'") = (1, . . . , 1) 
(to ones). It remains to recall that S(im) = Cm, the TOth elementary symmetric 
polynomial [Ml Ch. 1.3]. 

Let us perform the integration in the right-hand side of (4.7). We first integrate 
over z, and obtain: 



1 

(27ri)2 



dz 



dw 1 



N 

n 



W 

z - 



q"^ ' 



c{q"3-') 



c{oo) 



n 



1 



- a" 



1 

27ri J w 

c(oo) 



dw 



N+l-i 



Y[{w q'^ 



Next, the integral in w amounts to taking the coefficient by w ' in Jlr5itj(^ ~ 
qV-r-ry Recalling the generating series for the elementary symmetric polynomials 
[Ml Ch. 1.2] 

L L 

$^e™(2/i,...,yL)i™ = [](l + y,i), 



m=0 



r=l 



we see that the right-hand side of (4.7) is equal to that of (4.8). This concludes the 
proof. □ 
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we can compute 

the last ingredient used in the Eynard-Mehta type formula for the kernel qK (4.9). 
Namely, let us define for to = 0, 



4.5. Functions ipi{x \ m). Using the result of Proposition 4.10 

rd- 

,iV-l: 

I to) : = (<y9„i+i * ... * at * I N))(x) 
Proposition 4.11. We have for m — 0, . . . , N — 1: 



ijj^^x I m) 



1 



(27ri)2 



dz 



dw 



N 



W 



N+l-i 



-n 



c(oo) 

The contours in z and w are counter-clockwise and do not intersect. The contour 



€q{x) in z encircles the points , q 



.. (and not the points ^,q^ ^, 



.;. The 



contour c(oo) in w contains €q{x) and has a sufficiently large radius. 
Proof. Using Lemma [Xs] and the definition of I ^) ( |4.3| , we can write 

tpi{x I to) 



7n(y-x) iq^ '^^^g)Af-m-l 



y : y>x 



m(vj~j-x) 



iq;q)N-m-i 



3 ■ fi-J>a: 



My I N) 

-[,v(^)- 



q)N-m-l 

Using the double contour integral formula for the inverse Vandermonde matrix 
(Proposition 4.10), we readily see that the summation over j : Vj — j > x turns into 
the integration over z G €q{x). This concludes the proof. □ 

Remark 4.12. Observe that the double contour integral expression for ilji{x \ to) 
of Proposition |4.11| makes no sense for m = N. However, there is another similar 



expression for 'ipi{x \ N) which directly follows from (4.3) and Proposition 4.10 

w — q 



Mx I N) = 



1 



(27ri)2 



dw 



1 



n 



c(g^) c(oo) 



r—1 



Note that this expression does not use the global contour €.q{x), and here z belongs 
to a small contour around q^ instead. 

4.6. Getting formula for the kernel qK. The Eynard-Mehta type formalism (Bj 
Thm. 4.4] (see more references in { 4.2 1 produces a formula for the correlation kernel 
qK based on the representation of the measure qV]\r^i, as a product of determinants 
(Proposition 4.4). The formula looks as (see [4.3 and [4.5 for definitions of the 
objects involved): 

qK{xi,ni;x2,n2) = -v^^"^^"^) (x2, xi) (4.9) 

rii N 

+ ^^[G-*],j((^, * <^(^'"^))(mri,xi) • tPj{x2 I n2). 



The inverse transpose of the "Gram matrix" looks as (see Lemma 4.8): 
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Using this fact and Lemmas |4.5| and |4.6| we can write 

/ XI-X2 + 1. „\ 

qJ^[Xi,ni,X2,n2) — —in2<ni>^X2<xi<l ] ^ 

[Q\ 9jni-n2-l 

rii N 

+ E( n il-q'-'))q^'-'^'''i^N+l-^{x2\n2). 

i—l r—ni + 1 

Observe that the function 'ipN+i-i{x2 \ "-2) depends on i only via the term 
under the integral (Proposition 4.11 ). It follows that we can compress the sum over 
i into a (7-hypergeometric function as follows: 

711 N 



E( n (1-9'--')) (9^^^"') 



{q''-';q-')N-n, ■ 2Mq-\q'''-'\q''-^ I q-'-.q^'w-^). 



Plugging in the expression for ipN+i~i{x2 I "-2) (Proposition 4.111 without the 
factor which went inside 2^1, we get the following formula for the kernel: 

r X1-X2 + I. „\ 



qK{xi,ni]X2,n2) = ~ln2<n^lx2<x^q'""-^''^ ""^^ ' 



{q;q)rn-r 



iq^-'-:q-')N-n, 
(27ri)2 



dz 



dw q 



— n2X2 ^112 



x2</'i(g-\g"^ 



€q(x2) C(00) 



1 „ni-l. „N-1 I „-l 



W W — Z 



,xi -1^ i^q^ '■. g)jv-n2-i TT '^^q"'' 



iqiq)N-n2-i ^Ji z - q"- 



n 



Here it is essential that rt2 < — 1. For 712 — N, one could write another formula 
for the kernel qK{xi, ni; X2, ^2) based on Remark 4.12 but since the top row of the 
array X — {x™} (corresponding to m = N) is fixed, we do not need such a formula. 

Finally, changing the variables as w = wq~^^ and z = zq~^^ and renaming them 
back to w, z, we see that Theorem 14. 11 is established. 



5. Correlation kernel for uniformly random 
Gelfand-Tsetlin schemes 

5.1. The kernel. In this section we compute the correlation kernel for uniformly 
random Gelfand-Tsetlin schemes with arbitrary fixed top row v £ GTjy: 

Theorem 5.1. The correlation kernel K of the uniform measure PAr^jy on interlac- 
ing particle arrays {x™ : m = 1, . . . , A^, j = 1, . . . , m} with fixed top row = Vj — j 
(j — 1, . . . , N ) is given for 1 < ni < N , 1 < n2 < N — I, and xi, 2:2 G Z 6?/ 



K{xi,ni;X2,n2) = ~ln2<nAx2<xi 



1 



jxi - X2 + l)ni-K2-l 

(ni - ^2 - 1)! 

(z -X2 + l)N-n2 



{N -712 - 1)!' 



N 



/rv -NT ^ dz (h dw r TT 

(27rl)^ J J (W-Xi)N-ni + l w-z^^^ 

€{x2) c(oo) 



W -\- r — Vr 

Z A- r — Vr 



(5.1) 

The counter-clockwise contour ^{x2) in z contains the points a;2, 0:2 + 1, . • . , i^i — 1 
and not X2 — I, X2 — 2, . . . ,vi^ — N . The counter-clockwise contour c(oo) contains 
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£(x2) (without intersecting it) and has a radius sufficiently large to include all the 
points Xi, — 1, . . . , xi — {N — ni). 



By taking appropriate top row particles as in (2.3), we see that Theorem 5.1 
readily implies Theorem [l] 



the 9 t 1 limit in Theorem 



4.1 



5.2. Preliminaries. In the rest of this section we prove Theorem |5.1| by taking 

That is, we obtain K ( |5.1| ) as 

(5.2) 



K{xi,ni;x2,n2) = \iin„K{xi,ni;x2,n2). 

Note that since the measures qPn.i^ tend to Vn,^ as g f 1; the correlation functions 
of Pat,,/ must be limits of those of qPjv,i/- But because the correlation kernel of a 
point process is not defined uniquely]^ relation (5.2 1 for any two correlation kernels 



of] 



'jv.i/ and qFN,!^ does not a priori have to hold. 



Remark 5.2. We do not know if it is possible to establish (5.2 1 directly by looking 



at the asymptotics of the integrand in (4.1 1 and transforming the contours in some 
way. Our proof involves breaking the integral for qK in (4.1) into much smaller 
pieces and looking at their Taylor expansions at g = 1. In such expansions almost 
all terms disappear in the g t 1 limit. This is the reason why a rather complicated 
formula for the g-deformed kernel qK of Theorem 4.1 (with a g-hypergeometric 



function inside) turns into a simpler expression for K (|5.1|) 

Observe that the additional summand in qK simply tends to the corresponding 
summand in K: 

(g-i-^+i;g)„,_. 



in 



(a;i 



X2 



r 



"1— "2 — 1 



(g;g)ni-n2-i (ni-n2 + i)! 

(this a well-known property of the g-Pochhammer symbols) . The rest of this section 
is devoted to establishing the convergence of the remaining double contour integrals. 



5.3. Residues. The integrand in (4.1) has (possible) poles in the variable z at 
the points g''i~^~^i , . . . , g'^""^"^^ . For (5.2| to hold, the residue at every point 
z = q^i~i~^^ (j = 1, . . . , TV) must have a limit as g f 1- Indeed, such a residue can 
be isolated by taking linear combinations of qK{xi,ni;x2,n2) with various values 
of X2- Fix j = 1, . . . , N. The residue at z = g'^j~-?~^i looks as 



qReSj 



('?;'?) 



N-n2-l 



-1. 



(5.3) 



dw 

2tt\ J w 

c(oo) 



■w 



q'^r^r-x.y 



For shorter notation, set 
(ai,...,aAr_i) := {vi - 1 



Xi, 



. - j - xi, . . . ,i^N - N - xi) eZ 



N-l 



(5.4) 



■^For example, one can conjugate the kernel as K(xi,ni;X2,n2) 
with a nonvanishing function f{x,n), which does not affect the correlation functions i2.5l 



/(xi ,ni ) / \ 
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Now let US compute the integral in w in (5.3). This amounts to simply taking the 
free term in the product of the generating series 20i and nr=i^(''^ ~ q"''). Thus, 
our next goal is to understand the q 'f 1 asymptotics of the following expression: 



;g 

27ri 



dw 



-1 „N-1 



q ^■,w 



c(oo) 
rii-1 N 



^))(-l) 



CN-i-iiq"",- ■ 



(5.5) 



i—O r— ni + 1 

Denote the right-hand side of this expression by g£'„j (ai, . . . , apf^i). 

5.4. Taylor expansion of g-specialized elementary symmetric polynomi- 

1 in gEn^{ai,...,aN-i) ( [5^ to get its g f 1 

un-i. Thus, we 
at g = 1 to 

see which terms would matter in the q t 1 asymptotics in (5.5) and therefore in the 
residue qReSj (5.3 1. 

We will need to use the following symmetric polynomials: 



als. One cannot simply plug q 
limit; in particular, that would destroy the dependence on ai, . 
need to consider the whole Taylor expansion of eAr_i_i((7"i, . . . 



E 

«1,---,««(A) 



where the sum is taken over all pairwise distinct indices ii, . . . ,ig(^x) from 1 to 
— 1, and A G Y is a Young diagram [Ml Ch. I.I] with number of parts i{X). Here 
ylm _ 1) _ _ TO ^ 1) is tJie falling factorial power. These m^'s are the 

so-called augmented factorial monomial symmetric polynomials. 



Proposition 5.3. For all m = 0, . . . , N — 1, we have the following two expansions: 



,9 



|A|! 



q 



,ai + ...+Q!jv-i 



(g- 1)1^1 (N 



|A| 



m^(Q;i, . . . ,ajv-i) (5.6) 
m^(-ai, . . . , -aAr_i), (5.7) 



where the sums are taken over all Young diagrams A G Y 
rii 

r,aN-l ^ 



To see that both series in ( 5.6 ) and converge, one could argue as follows. Us- 



ing em.{q° 



one could make all a^'s positive. Since m^(/3i. 



(cm's are homogeneous), 
Pn-i) vanishes for nonnegative 



/3i's if all Aj's are large enough, the sum in (5.6) for positive a^'s will be finite. 
Thus, one can turn (5.6) into a power of q (which expands into a convergent series 
at g = 1) times a polynomial in g — 1. The same trick can be performed for (5.7). 



Proof. The two expansions (5.6) and (5.7) turn one into another because 
,yN-i) = yi ■ --yN-i ■ eJv-l-m(2/^^ • ■ 



i(yi 



-yAT-i)- 



Thus, we only need to prove, say, (5.6). 
Consider the generating function 

N-l 

Fit;q) := ^i"e™(g"S. 

m— 



N-l 



■ ,q 



Qn-1 \ 



)= 11(1 + *?"'■)• 
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Since 

^(l + ^(Z"'■)|,=l=toi:^ s>l, 
we have for any s = 0, 1, . . . the following expression for partial derivatives of F{t; q): 

^F{t; = ^'^'^(1 + tf-'-'^'^^lic^i, c^N-i). (5.8) 

Indeed, every summand corresponding to some partition A indicates that we ap- 
ply the derivative to one of the factors (1 + tq"''i) for every i = 1, . . . ,i{X). 
Summing over all possibilities of doing that, we get the polynomial m^. After 
this, — 1 — ^(A) factors are not affected by differentiatio n, an d this gives us the 
multiplication by {l + t)^~^^~^^^K The special case s = in (5.81 is checked directly. 
Since 

iV - 1 - £(A) 



(coefficient by t™), we obtain ( |5.6| by writing the standard Taylor expansion. This 
concludes the proof. □ 

5.5. g-Stirling and classical Stirling numbers. Now, plugging the expansion 
of eN-i-i{q"^ , . . . , q""-^) given by (|5.7| into (|5.5[), we can rewrite 



X 

i—O r— ni+1 



In this subsection we will relate the sum over i above to the known g-Stirling 
numbers (of the second kind) , and this will help us to write the leading term of the 
q t 1 asymptotics of g-E„j (ai, . . . , ctN-i) (Proposition 5.7 below). 
We will use the standard g-notation: 

1 - (j" 

[n]q := ^ N-?' ■= W'?!'^ " 1]« ■ ■ • W?- 

Definition 5.4 ( [Gouj ) . The q-Stirling numbers of the second kind S{n,m;q) are 
defined as the following expansion coefficients (here l^l is sufficiently small): 

n oo 

n i_H ,^ -T.S{n,m;q)z-^. (5.10) 

r=l ^ -I"? m=0 

Lemma 5.5. For every I ~ 0, . . . , N — 1, we have 

Jil — 1 N / AT T I\ 

E( n (s-^d 

■i=0 r=ni-|-l ^ 

= (-l)^+«i (1 - q)^-'-^[N - m],! • S{N - n,,n, - I - 1; q). 
Proof. We use the following explicit formula for the q-Stirling numbers from |Gou| : 



5(n,m;g) = (Q-l)-"^(-l)f 

p=0 



plm + n\ [m + n~ p]ql 



]gl[m-p]gl' 
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Thus, 



S{N — rti, rti — / — 1; 5) 

— ^ n 1 — / — 1 



= (g-i 



(9 



Til — 1 



V 

N-l- 1 



[N-1-1-P]g\ 



[N - ni]gl[ni - I - 1 - p]gl 
[TV-l-p],! 



7V-p-iy [7V-ni]g![ni-l-p], 



Note that the summation over i in (5.11) is also done from Z to ni — 1 due to the 
presence of the binomial coefficient (^Z}z')- It remains to match terms in the 
above formula with those in ( 5.11[ ) to see that the claim holds. □ 

The g-Stirling numbers S{n,m;q) converge to the classical Stirling numbers 
S{n, m) which are defined as the following expansion coefficients (for sufhciently 
small \z\): 



n 



00 

E 

m=0 



S{n, m)z'" 



(5.12) 



Remark 5.6. It can be readily checked that the numbers S{n, m) are related to 
the more common parametrization of the Stirling numbers of the second kind 



1 



ml ^ 



4f 



(to - p)" 



via S{n,m) = { "^'"} • 

Proposition 5.7. The leading term of the asymptotics of qEn^ (ai, . . . ^a^- 
as q ^ 1 looks as follows: 

,i?„,(ai, . . . , aN-i) ~ (1 - qf-^-lf+'^HN 



ni y.x 



, aN-i)S{N — ni,ni — I 



1) g 

(5.13) 
!)• 



1=0 



Pro of. Looking at formula (5.91 and Lemma 5.5 we see that each sum over i 
in (5.9) behaves as ~ const • (1 — g)^^^'"^^^^. Therefore, the whole expression 
qEn^ {ai, . . . , Q!jv_i) (5.9 ) behaves as 



const • (1 — g) 



constA- (<z-l)l^l"^^^^ 



It follows that Young diagrams A with |A| > £{X) provide a negligible contribution 
to qi5„j(ai, . . . ,ajv-i). Thus, we arc left only with one-column Young diagrams, 
and for each of them the factorial monomial symmetric polynomial reduces to a 
multiple of££nv) . Lemma 5.5 then provides the necessary asymptotics which directly 
leads to (5.131. Observe that we only need to sum over < Z < ni — 1 because 
for bigger /, the Stirling numbers S{N — ni,ni — I — 1) vanish. This concludes the 
proof. □ 



Remark 5.8. We see from the above proposition that taking g t 1 limit allows us 
to drop almost all summands in the sum over A € Y in (5.9 1. This is the reason why 
the g-hypergeometric function in the formula for the kernel qK (4.1 ) disappears for 
q~l (5.1). Proposition 5.9 below also displays this effect. 
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5.6. Completing the proof. Summarizing the development of § §5.4f|5.5| and trans- 
lating Proposition |5 . 7| into the language of contour integrals, we have: 



Proposition 5.9. The leading term as q'\ 1 of the contour integral in w in qReSj 

dw 



{5.3) looks as follows: 

-N-l. -1 



c(oo) 



q ^■,w' 



'')l[iw-q 



Vr — r—xi 



(g-l)^-^(iV-ni)! 



dw 



c(oo) {W - Xi)N-m 



+1 



n( 



w - 



(5.14) 



Here in both integrals the contours in w are counter-clockwise and have sufficiently 
large radii. 



Proof. A simple change of variables in the integral in the right-hand side of (5.14) 
and the expansion of the integrand into powers of w lead to appearance of the 
classical Stirling numbers via (5.121, and also of the elementary symmetric polyno- 
mials ei{ai, . . . , aAr_i) (recall (5.4)). Then it can be readily checked that the claim 
directly follows from (5.5) and Proposition 5.7 □ 



Now we are in a position to compute the limits of the residues qResj (5.3) 

(iV-ni)! {l^j - j - X2 + l)N-n2-l 
1)! Ur^M 



lim „i?es, —: Res^ — 



{N - 712 



27ri 



c(oo) 



dw 

{w - Xi)Ar_„i+i 



X (5.15) 



5.9 



Indeed, the factor (q — 1)^ ^ in Proposition 
the product Hr^ji^l''' 

Everything else in the above formula also follows from (5.3) and Proposition |5. 9 



3-^^ - q^^-'-^^)-^ in"]5.3[) into n.^j( 



is exactly what is needed to turn 



To complete the proof of Theorcm |5.1[ it remains to note that the double contour 
integral in the formula for the q-deformed kernel qK{xi,ni; X2, f^2) (4.1 ) is equal to 
the sum of qReSj ( 5.3 ) over all j = 1 , . . . , TV such that Vj — j > X2', and the same is 
true for the double contour integral for K{xi,ni; X2,n2) in (5.1) and the limiting 
residues Resj defined above. 

Thus, we have completed the proof of Theorem |5.1[ and also of Theorem [T] 



6. Inverse Kasteleyn matrix 



6.1. Inverse Kasteleyn matrix and the correlation kernel. Let us show how 
the kernel K is related to the inverse of the Kasteleyn matrix for the honeycomb 
graph Gp inside our polygon P (Fig. [l] right). We would like to use the affine 
transform (2.1 ), so that P will be parametrized by Ai, Bi as in i|2l](see Fig.|3|, or, 
equivalently (see (2.3) and ^3.5), by a fixed signature v e GTn as in Theorem |5.1| 
The graph Gp is bipartite; its vertices correspond to two types of (triangle) faces 
in the dual triangular lattice: 
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We will encode each such triangle by the position (x, n) of the mid-point of its 
horizontal side. The Kasteleyn matrix of the graph Gp is its adjacency matrix 
with rows and columns parametrized by white and black triangles, respectively 
(e.g., see |K3| ) . Inside the polygon, this matrix looks as 



Kast(^(x, n); L.(y, rn)) 



(6.1) 



'1, if (y,m) = (a;,n); 
1, if (y,m) = (a;,n- 1); 
1, if (y, m) = (x + l,n — 1); 
0, otherwise 

(see Fig. [9]). For ^(x, n) on the boundary of the graph Gp, the ^{x, n)-th row of 



Figure 9. Edges of three directions in the graph Gp encoded by 
pairs of triangles. 

Kast will contain less than three ones, and the same for the ^{y, m)-th column. 

The inverse matrix Kast~ has rows and columns indexed by black and white 
triangles, respectively. As |K3[ Cor. 3] suggests, Kast~^ can also serve as a correla- 
tion kern el fo r the uniform measure on tilings of P. Based on the explicit formula of 
Theorem 5.1 (or Theorem[l]), we establish the following connection between Kast^^ 
and our correlation kernel K: 

Theorem 6.1. The inverse Kasteleyn matrix and the correlation kernel K of The- 
orem \571\ are related as follows: 

Kast-i(k(y,m);^(a;,n)) = (-l)^-^+™-"iC(a;, n; y, m). 
Proof. We need to check that 



(-l)^'""+™'-"i^(x,n;y',m')Kast(^(a;,n);k(y,m)) = l(,,„)=(^-,„o- (6-2) 



The black triangles L.(?/, m) in the graph Gp have < m < — 1 (e.g., see 
Fig.|3|. Depending on the position of a black triangle, (6.2 ) turns into a three-term 
or a two-term relation for the correlation kernel K . We aim to verify all these 



relations using the explicit formula (5.11 which expresses the correlation kernel 



-fir(a;i, rii; X2, ^2) as a double contour integral (denote it by /(xi, ni; a;2, 7^2)) plus 
an additional summand. 

1. (inside the polygon) Let 1 < rn < A — 1, and y be such that the L.(y, m)-th 
column of Kast has three ones. Then \<o.2\ becomes the three-term relation: 



K{y,m;y\m') ~ K{y,m + l;?/',m') + K{y - l,m+ l;y',m') = l{y,m)={y' ,m')- 

(6.3) 

Using an obvious identity 

{N~my. {N-m-iy. {N - m - 1)1 



{w - y)N~ 



-m+l 



(w - y)N- 



{w -y+ l)N-ni 



0, 
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we see that the terms corresponding to the double contour integral part of K sum 
to zero: 

/(y, m; y\ m') ~ I{y, m + 1; y' , m') + I{y-l,m + 1; y' , m') = 0, (6.4) 



because the contours in all these integrals are the same. Note that (6.4) in fact 
holds for any < m, m' < TV — 1 and any y, y' £ Z. 

The fact that the terms corresponding to the additional summand in K give 
the desired result l(j,,m)=(i/',m') in (6.3) can be checked directly. Thus, (6.2) is 
established inside the polygon. 

2. (boundary) Now we consider boundary black triangles ^.(y, m). Note that the 
ones which are close to the top boundary of P are already included in the general 
case 1. 

2a. (bottom boundary) Let to = 0, so the triangle L.(y, to) lies at the bottom 
of P (Fig. [3]). For to = 0, turns into 

-K{y, 1; y', m') + K{y - 1, 1; y', to') = l(y,o)=(y',m')- (6-5) 

One can show that _ftr(y, 0; y', to') vanishes: (1) there is no additional summand in 
(5.1) because m' > 0; (2) looking at the integral in w G c(oo) in (5.1), we see that 
the w-integrand decays as at oo and thus has no residue there, so the integral 



over w is also zero. Thus, we may include K(y, 0; y', to') into ( |6.5[ ), and use the fact 
that the resulting three-term relation (6.3) is already established. 

2b. (vertical boundary) Let 1 < to < — 1, and let L.(y,TO) be close to the 
vertical boundary of the polygon P (except for the rightmost vertical boundary 
which falls into the general case 1). This means that the point (y — 1, to + 1) is 
outside P, and we must show that 

K{y,m;y',m') ~ K{y,m + l;y',m') = l{y^m)=(y' .m')- (6.6) 



Using (6.4|, we can replace the double contour integral parts I{y,m;y' ,m') — 
I{y, TO + 1; y', to') above by — /(y — 1, to + 1; y', m'). 

We can compute the resulting double contour integral I{y — l,m + l;y',m'). 
Observe that for the point (y — l,m + 1) to be outside P, the top row particles 
must occupy positions y ~ N + m, . . . ,y — 1 (see Fig. [s]). This means that the 
integrand in w in I{y — l,m + l;y',TO') does not have poles inside the contour 
except for w ^ z. Taking the residue at w — z and using the Lemma 6.2 below, we 
see that 



I{y - 1, TO + 1; y', to') = ly,^y 



[m — m' + l)y-y' 



Taking into account the additional summands that come from the two kernels in 
(6.6), one can directly check that the desired identity (6.61 holds for every {y',m') 
inside the polygon P. In fact, (6.6) may fail if (y', m') is outside the polygon. 

2c. (boundary parallel to the vector (—1, 1)) Finally, let 1 < to < — 1, and 
L.(y,TO) be close to the boundary of P of direction (—1, 1) (except for the leftmost 
such boundary which falls into the general case 1). We need to show that 

K{y, to; y', to') + K{y - 1, to + 1; y', m') = l(y.m)=(y' ,7n') ■ 

Since L.(y, to) is at the boundary, and thus the point {y,m + 1) is outside P, we 
note that the top row particles must occupy positions y — N -\- m — l,...,y (see 
Fig.[3|. Then one can argue in the same way as in 2b, and also use Lemma |6.2| 



This concludes the proof of the theorem modulo the below Lemma p?2\ □ 
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6.2. Computing some contour integrals. 

Lemma 6.2. For all 1 < ni < N , I < n2 < N — 1, and xi,X2 € Z, we have 



{N-n2-iy. 



{z- X2 + l)Ar-„2-l 



. (ni - n2)xi-x2 If. 

{x,^x2y. ' ^^-^^ 



where the contour €,{x2) in z is the same as in Theorem 5.1' a counter-clockwise 
contour which encircles points X2,X2 -\- 1, . . ■ , and not points X2 — 1^X2 — 2, . . .. 



As was mentioned in tlic proof of tlie above Theorem 6.1 tlie integral in (6.7) 
arises if one takes the w = z residue in the w integral in the formula for our 
correlation kernel K{xi,ni] X2,n2) (5.11. 



Proof. Denote the integral in (6.7) by J. Let us denote Ax = xi — X2 and An 
ni — 7i2. Shifting the variable z hy Xi, we have 

f {z + Ax +l)N-n2-l 



J = 



{N-ni)\ 1 
(TV - 712 - 1)! ^ 2^ 



-dz. 



{z)N-n-i_ + l 

The integrand here has poles 0, —1, . . . , —{N — ni). We see that if Aa; < 0, the 
contour of integration €.{—Ax) has no poles inside, and thus J = 0. 

In the rest of the proof we assume that Ax > 0. The integral J then becomes 
a sum over the poles 0, —1, . . . , —Ax of the corresponding residues which we write 
as follows: 

r(Aa; + l) r(A^-ni + l) 



J : 



r(iV - na + Ax) 
T{Ax + 1)T{N - n2 



E 

J=0 



r{z + Ax + 1) r(z + - ni + 1) 
r{z + Ax + N - n2 



V{Ax + N -n2) 



{z + j)T{z) 



Using a simple observation that 

T{-A+l) T{-A+l+j-j) 



T{-A+\-]) 



T{-A + 



{-A+l-j),^{-iy{A), 



(6.8) 



and the residue of the Gamma function (2+ j)r(z)|z=_j — (— l)^/j!, we can rewrite 
as follows: 



J = 



r(7V - + Ax) 



Ax 

E 



{-Ax),{ni- N)j 



1 



r(A.T + l)V{N - n2) ^ {-Ax -N + n2 + I), f 
r{N -n2 + Ax) 



2F1 



-Ax,ni - N 
-Ax - N + n2 + l 



T{Ax + l)r{N-n2)' 
Here 2^"! is the Gauss hypergeometric function. We may use the Gauss summation 
formula [El 2.8.(46)] for it if we assume for a while that iV is a nonreal complex 
number: 



-Ax,ni - N 

-Ax - N + 712 + 1 



T{-Ax - N + n2 + l) r(-An + l) 
r{-N + n2 + 1) r(-Aa; - An + 1) 

(An)Aa: 

{N - n2)Ax 



(we have used (6.8) again). Then we can set the complex N to be equal to an 
integer again because both sides of the above identity are rational functions in N. 
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Putting all together, we have 

r{N-n2 + Ax) (An) Ax 



J 



T{Ax + l)T{N-n2){N-n2)Ax r(Aa; + 1) 



This concludes the proof of the lemma and completes the proof of Theorem 6.1 □ 

Let us compute several related contour integrals which will be useful for the 
asymptotics at the edge (Sjsl). 



(6.9) 



Lemma 6.3. For all 1 < rii < N , 1 < n2 < N ~ 1, and Xi,X2 € we have 

{N-my. ^ J_ / {Z-X2 + l)jv-»2-i ^^ 
(iV-n2-l)! 27ri/c,(^2) {z - xij^^rn+i 



>x- 



- 1 



ni>n2/ 



ni — n2 + Xi — X2 — 1 

Xi - X2 



where the contour £'(2:2) in z is a clockwise contour which encircles points X2 
1, X2 — 2, . . and not points X2,X2 + 1, ■ ■ ■■ 



Proof. The proof is analogous to that of Lemma 6.2 



□ 



Note that for the two integrals in (6.7) and (6.9) we have /g^^^^ — .^£'(2:2) ^ .^c(oo)' 
where the last integral is taken over a large enough counter-clockwise contour. Since 
in (6.7) and (6.9) the binomial coefhcients in the right-hand sides are the same, we 

{z-X2 + l)jV-«2-l 



also have 
{N-n,y. 



1 



(TV -77,2-1)! 27ri / {z - Xi)N-ni+l 
c(oo) 



-dz = 1 



{xi — X2 + l)ni-n2-l 



rii >n2 " 



(ni - 712 - 1)! 



(6.10) 



Lemma 6.4. For all 1 < ni < N , I < n2 < N — 1, and xi,X2 G we have 

N 



1 



{N-n2-iy (27ri)2 



dz 



dw 



{z- X2 + l)N-n2-l 1 



n 



w + r — l^r 



c(oo) c(oo) 



{W - Xi)N-ni + l ^ Z + r-l^r 

_ -I (xi -X2 + l)ni-n2-l 

-l"^>"^ (m- 772-1)! ■ ^^-^^^ 



This is the same integral as in the correlation kernel {5.1), but with the contour 
£(2:2) in z replaced by a sufficiently large counter-clockwise contour which lies inside 
the w contour. 



Proof. Let us write the integral in z in (6.11 ) as a sum over the residues a.t z = Vj—j, 
j = l,...,N: 

1 



(27ri)2 



dz 



dw 



[z - X2 + l)N-n2-l 1 



N 

n 



c(oo) c(oo) 



27ri 



w + r — Vj- 

{w-Xi)N-m+i w - z ^J-^ z + r - Vr 

N 



w + r — Vj. 



— / 7 — — y'(i^j-j-2;2+i)Ar-n2-in 

/ (Ti; - a;i)Ar_„^+i ^ l^r, ^3 - J + r - i^r 

c(oo) ■' 
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The sum here is the Lagrange interpolation polynomial of the function w t-^ (w — 
3^2 + l)Ar-n2-i with nodes — — 1,...,N. Since that function in w is itself a 
polynomial of degree < iV — 1, the interpolation is exact, and 

Wj - 3 - X2 + IjAT-na-l — ^ [W - X2 + l)w-«2-l- 

j = l r^j J J ^ 



The claim now follows from (6.101. □ 



7. ASYMPTOTICS IN THE BULK, LIMIT SHAPE, AND FROZEN BOUNDARY 

7.1. Parameters of the polygon. In this and the next section we perform asymp- 
totic analysis of the uniform measure Pp(7v) on lozenge tilings of the polygon P(Af) 
as iV oo. We will use the parameters k and {Ai{N),Bi{N)}\^-^ C Z' = Z + i 
describing the polygon P(A^) ([p]). We will think that k = 2,3,... is fixed, and 



Ai{N), Bi{N) scale linearly with N as in (2.7), depending on new continuous pa- 
rameters {ai,?>i}i=i with oi < 6i < . . . < afe < &fe, and ~ ^i) — 1- These 
parameters describe the limiting polygon V in the new coordinates {x,il) (Fig. |4|. 

We fix a global position (x, rj) e V, and obtain local asymptotics of the measures 
Pp(jv) around (x, ^7) via asymptotic analysis of the correlation kernel K of Theorem 
[1] We then discuss how these local asymptotics describe global properties of random 
lozenge tilings such as the limit shape and the frozen boundary. In this section we 
prove Theorem [2] and Propositions [231 [2]6) and[2J] 

7.2. Asymptotic expression for the kernel in the bulk regime. As a first 
step, we establish an asymptotically equivalent expression for the kernel K which 
will allow to employ saddle point analysis in the spirit of |Okj (see also [ORlj . 
[QR2], [BK]). 

We will look at the kernel K{xi,ni; X2,n2) with the parameters scaled as 

xi^2/N^X, ni^2/N^r], A ^ 00, (7.1) 

such that the differences 

Ax :— xi ~ X2 (z Z, An := ni — 712 G Z (7-2) 

stabilize. This is the so-called 'bulk' limit regime which describes local asymptotic 
behavior of the measures Pp(Ar) around the global position (XtV)- 

Definition 7.1. Define the action by 

S{w; X, v) ■= {w - x) ln(w - x) - (w - x + 1 - »y) ln(w - X + 1 - ?/) (7-3) 

k 

+ {l-Tj) ln(l - f?) + ^ [(foj - w) ln(6i - w) - {ai - w) ln(ai - w) . 

i=l 

Unless otherwise stated, we assume that that the branches of all logarithms have 
cuts looking in negative direction along the real line. Note that the real part 
'^S{w;Xt1) is well-defined and continuous for all w G C. 



Proposition 7.2. In the regime (7.1)-(7.2), the correlation kernel K of Theorem^ 
is asymptotically equivalent to 

u( ^ 1 1 (Aa: + l)An-i ^ (l->?)^'^" 1,1, 
K{xi,ni;x2,n2) ^ -lAn>o^Ax>o — _ 1 ^^^.-^2 f f 

ff(x-) c(oo) 
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X ^ — ^ . (7.4) 

(z-x)5(^-X + l-J?)5 w ~ z 

The branches of the square roots and other noninteger powers here are assumed to 
have cuts looking in negative direction along the real line. 



The z contour €{x~) i''^ l''^-4'l counter-clockwise, it starts inside the segment 
(x + ?7 — goes in the upper half plane, crosses the real line again to the right 

of bkH and returns (in the lower half plane) back to where it started. The counter- 
clockwise w contour c(oo) contains ) without intersecting it, and has a suffi- 
ciently large radius. See contours on Figure [771 ^^ft- 



Note that by our choice of branches, the integrand in (7.4) is continuous on our 
contours except for possibly one real point lying on a cut. 
In the rest of this subsection we prove Proposition |7.2[ 



By scaling the variables of integration as z — z/N, w = w/N in formula (2.6) 
for K{xi, Hi] X2, (and renaming back), we can write 

irf \ -I ^ (Aa; + l)A„-i 

K{Xi,ni;X2,n2) = -lAn>olA2;>0 7T TTj — (7.5) 

^ (An — 1)! 



(2^i)2 J J + P(^; 3:2,^2)' 

ff(X-) c(oo) 

where 

P{w;x,n) := constAr • — — — '■ TT + 5 ^ Nw)B,~Ai- 

(Nw -X-\- l)N-n-l fj^ 

Here and below constAr denotes a constant which does not depend on w,x,n, but 
may depend on A^. Clearly, such a constant in P does not affect the integrand in 



(7.5) 



Let us explain why we can choose the contour £(%— ) for z in (7.5 1. After the 



scaling, the new z contour will contain inside it the points =^ , , . . . ^'''^^^ , and 
the points ■ • ■ will be outside. It is possible to drag the left end of the 

contour slightly to the left because the factor {Nz — X2 + 1) Ar-n2-i inside 



P(z;X2 ,"2) 

in (|7.5|) compensates the corresponding poles. We can also drag the right end of the 



contour slightly to the right. Thus, we arrive at the macroscopic contour £(%— ) for 



z (it does not depend on N). The z and w contours in (7.5 1 are given on Fig. 11 
left. 

Lemma 7.3. Let x/N — > x and n/N — 77 as N 00. Then 
P{w- X, n) = const^ ( ~ ^^^"^ ~ | ^ ^ " ) ' exp (nS{w- ^, ^) + 0(i)) 



for w G C \ M. Here S is given by {7.3) 



Proof. We write every Pochhammer symbol in the definition of P(w; x, n) as a ratio 
of two Gamma functions, and then apply the Stirling approximation (which holds 
for y ^ (-00, 0]): 

r(y)-exp((y-i)ln2/-y+iln(2^) + 0(i)), |y| ^ 00. (7.6) 



"^Note that (x, »?) G P in particular means that X < ''fc ■ 
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The desired asymptotics then fohows from a direct computation involving Taylor 
expansions of logarithms arising. □ 

Let us now deal with the exponential factor in -^('"■^i'"!) ■ 

^ F(z;X2,n2) 



Lemma 7.4. In the regime (7.1)-[7.2), we have for w, z ^ C\ 



ri2 'lAx+Aw 
N ) 



exp 



(n[S{w- 



N ^ N I N ^ N l\ )■ 



Proof. Indeed, we can write 



Ax 712 

N ' N 



An 
N 



_ Q(oir ^2 "2 \ I Ax c X2 n2 ^) I An c („,,. X2 112 \ 

where the derivatives are given by 

Sx(w, X, v) = w^S^iw; X, 77) = ln(w - x + 1 - '7) - - x), 
S,,{w; X, v) = |^'5'x('^; X, V) ^ ln(w - X + 1 - - ln(l - ??)• 

This concludes the proof. 



□ 



7.3 



To see that Proposition 7.2 holds, it remains to put together (|7.5|) and Lemmas 
and 



7.4 



Outside the exponent, we clearly can replace and by x and 
77, respectively. The fact that uijZ G C \ K in Lemmas 7.3 and 7.4 affects only two 
points on each of the contours, and changing values of the integrand at these points 
does not change the integral. Thus, Proposition |7.2| is established. 

7.3. Critical points of the action S{w,x,v)- As the saddle point technique 
suggests [Ok] , to analyse the asymptotics of the double contour integral for the pre- 
limit kernel (7.4), we need to deal with the critical points of the action S{w;x^v) 
(Definition 7.1), i.e., solutions to 



-,S{w,x,v) = 0. 



(7.7) 



Proposition 7.5. For every point (XtV) inside the polygon V (defined by {0^,6^}, 
see Fig. the function S{w, Xj ""^ ^'^•s either 2 or nonreal critical points. 



Proof. The equation ( |7.7[ ) for the critical points is equivalent to the following alge- 
braic equation of degree k: 

{w-x)X^. iw-a^) = {'w-x+l-ri)\^, iw-b,). (7.8) 

Let us denote by Qa{w) and Qb{w) the polynomials on the left and on the right, 
respectively. We aim to count the number of intersections of graphs of Qa{w) and 
Qb{w), w e M. Clearly, 

Qa{w) — Qb{w) = rjw^^^ + lower degree terms. 

Inside the polygon V , one has < 77 < 1, so we know which of the two graphs of 
Qa{w) and Qb(w) is asymptotically higher dX w — ±00. Thus, knowing roots of 
Qa{w) and Qb{w), we can directly count the number of intersections (an example 
shown on Fig. 10), and see that the proposition holds. □ 



ASYMPTOTICS OF RANDOM LOZENGE TILINGS VIA GELFAND-TSETLIN SCHEMES 33 




Figure 10. Example with k — A. The polynomials Qa{w) and 
Qb{w) intersect 4 (left; the fourth intersection point is far to the 
right) or 2 (right) times depending on (x^v)- 



Remark 7.6. In figures below where we show contours of integration or other 
relevant objects, the segments [oi, bi] and the distinguished segment [x + V ~ 1; x] 
are also present. We will always follow the way these segments are shown on Fig. \T0\ 
there are several red segments [ai,bi] lying on the horizontal line, and one black 
segment [x + ^7 ~ liX] which is shown slightly under that line. Its the endpoints 
X + — 1 and x are indicated by small vertical marks crossing the horizontal line. 



Definition 7.7. According to Proposition 7.5 \ct V d V be the (in fact, open) 
set of pairs (XtV) for which S(w;XjV) has nonreal critical points. This set V is 
called the liquid region, it lies inside the frozen boundary curve dD. See Fig. [4] 
and § plH2^ 

For {XjV) G let Wc = Wc{x,il) denote the unique nonreal critical point of 
S{'w; Xj v) lying in the upper half plane. 

7.4. Moving the contours. Fix a point {x, rj) G T). We aim to move the contours 



in the double contour integral for the correlation kernel if (xi, ni; X2, ^2) (7.4) so 
that the exponent exp |^iV[S'(w; j^, j^) — S{z; 7/)]) will make the whole inte- 
gral go to zero. This is achieved by making the z and w contours cross at the two 
simple nonreal critical points Wc, Wc of S{w; Xj v)- These points are also the saddle 
points of 5ft5(w; x, v)j so on the new transformed contours we have 

^S{w;j^,^)<nS{w,-x,v)<^S{z;^,^), z,w^^n, (7.9) 

(see, e.g., |0k[ §3] for more detail). In the course of moving the contours, certain 
residues corresponding to w = z will appear. 

Proposition 7.8. It is always possible to transform the w and z contours in the 



double contour integral for K(xi,ni] X2,n2) {7-4) such that {7.9) holds (Fig. 11 
right). After this transformation, an additional term that converges to 

i^J^ / [z - xr^^^-^z - X + 1 - r;)^-+^"-idz (7.10) 



will arise. In (7.10), the contour of integration crosses (x, +00). 



Proof. It is clear from the original formula for the kernel K (2.6) that one can drag 
the z contour C(x— ) everywhere along the real line except for the segments [oi, bi] 
which are not covered by [x + ^7 ^ I7 x]- Likewise, one can drag the w contour c(oo) 
everywhere except for the part of [x + — 1, x] outside the segments [oj, bi]. 

Observe that for any (x, rf) G V, the real part diS{w; x, v) goes to +00 as |w| — ^ 00 
because ^S{'w; x, r]) ^ rjln ]w]. This implies that the only possible configuration of 
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Figure 11. The contours before (left) and after (right) transfor- 
mation. Shaded is the region where ^{S{w; x, v) ~ S{wc', XiV)) > 0- 



shaded regions around the saddle points is as on Fig.jll] This means that we can 
let the z contour spread to the right arbitrarily far. 

Dragging the w contour through the z one, in the w integral in ( 7.4 ) we pick up 
the residue at w — z for z between Wc and Wc. Then we integrate this residue over 
z, which results in (7.101. 

It remains to show that there exist: 

• a point zi inside [x + V ~ IjX] where 3?S'(w;x, J?) > ^S{\Nc; x^v)- This is the 
desired left point of intersection of the z contour with the real line; 

• two points wi and Wr outside [x + V ^ IjX] where the opposite inequality holds. 
These are the desired points of intersection of the w contour and the real line. 

Existence of such points can be established by looking at the imaginary parts 
'^S{w;XjV) ^-iid 5S'(wc;Xj?7) and finding the points w where they are equal. We 
assume that now the cuts in the definition of S (7.3) are in the lower half plane, 
so 5(S'(w;x5??) — S{wc]XjV)) is continuous in the upper half plane up to the real 
line (except for a finite number of points). Counting the arguments, one sees that 



Figure 12. Curves along which 5(S'(w;x, 77) - S{\Nc;XiV)) = 0- 
Each such curve from to Wj, curve completely belongs to a 
shaded or not shaded region because of the local behavior of 
S{w;XtV) around its critical points. 

there is one point inside [x + '/ ^ 1; x]i and two points outside that segment where 
'^■{S{w; X, v)~S{\Nc; XiV)) = 0- Wc can choose them to be the desired zi and wi,Wr, 
respectively, see Fig. [12] □ 

7.5. Proof of Theorem[2} Now we can finalize the proof of Theorem[2]formulated 



in [2.3 Recall the incomplete beta kernel BQ(m, /) with complex parameter G 
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C, ^fi > (Definition 2.4). Wc will prove the convergence of the correlation 
kernel K{xi,ni] X2,n2) (2.61 to the incomplete beta kernel which will readily imply 
Theorem [2] via ( p^ . 



Proposition 7.9. For {x,v) ^ *^ bulk limit regime (7.1)-(7.2), we have the 
following convergence: 

K{xi,ni;x2,n2) 60(712 - ni,xi - X2), 

where the parameter = Q{xi'n) 0/ the incomplete beta kernel is given by 



(7.11) 



Wc(x,'7)-X+1-'?' 
Proof. Transforming the contours as in Proposition |7.8[ with the help of Proposition 



7.21 we see that 

K{xi,ni;x2,n2) 



-1 



An>olAa;>0 



(1-77)1 



-An 



27ri 



(Aa; + l)A«-i 
(Ati- 1)! 

(z - x)-^"- 



-dz 



(27ri)2 



(z - X + 1 - 77)-^^-^"+! 
(w - x)"'^"'"^ (w - X + 1 - ?y)'^^+^«-i 



(2^ - x) - X + 1 - ??) ^ 



-SO 



^ "2 X 



ui — z 

The first two summands give Bn(— An, Aa;): (1) for An > 0, drag the contour in 

A^^+^)^^-^ which wiU 



/- " through the point x, this will result in the residue 1ax>o (An-^y. 
compensate the first summand; (2) in the integral substitute u — ^ J^^^_.^ , this 
will give the integral for Bq of Definition |2.4[ 

The contours in the third summand above are as on Fig. |11[ right, and so it goes 
to zero because of the exponential factor]^ This concludes the proof. □ 

For (XtV) i in the bulk regime, the kernel K{xi,ni; X2,n2) also converges to 
the incomplete beta kernel B^^, but with real fl. This can be seen also by moving 
the contours as in Proposition |7.8[ but now all saddle points would be located 
on the real line. One of the transformed contours should pass through a saddle 
point, and it will contain or will be contained inside the other contour. The residue 



picked up after the transformation will be either 0, or the integral (7.10) over a 



closed contour. This exactly corresponds to with real il. A detailed treatment 
of contours for the case of real saddle points for a similar double contour integral 
is performed in |BK1 §4.4]. 

We have now established Theorem [2l 

7.6. Limit shape and the complex Burgers equation. Here we prove Proposi- 

it readily follows that the complex slope = fl{x, v) 



tion 



2.5 From (7.111 and (7 



satisfies the algebraic equation 

k 

n-Y[{{a,-x + i-v)^- 



^One should take the saddle points Wc(x2/N,n2/N) for each A^th integral, and transform the 
contours accordingly. But since X2/N x, n2/N rj, these saddle points are close to Wc(x,7?)- 
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This equation has degree fc + 1 in contrast with the equation ( 7.8 1 for Wc which has 



degree fc. However, the above equation for Q always has a root i7 = 1, so it can be 
reduced to a certain degree fc equation. 

Proposition 7.10 (Complex Burgers equation). In addition to the above algebraic 
equation, il{x,'r]) '^^■so satisfies the complex Burgers equation |KOj ; 

^{x.m — — ^ -(1 - ^{x.m) — — • 

Proof. It is more suitable to start with a differential equation for the parame- 



ter Wc(x,'7) of Definition 7.7 The equation for the critical points of the action 



S(w\ X, 77) ( |7.3[ ) which is satisfied by Wc can be written as 

fc 

ln(w — x) ~ ln(u' — X + 1 ~ ^) = ( ln(foi — w) — ln(ai — w)) 



(7.12) 



1=1 



(it is of course equivalent to (7.8)). Differentiating it with respect to x and 77, we 
obtain: 



1 



(Wc)> 



1 



Wc - X 
(Wc)r, 



Wc - X + 1 - 



(Wc)x ■ S(Wc), 



(Wc), 



1 



Wc-X Wc-x + l- ?7 



= (Wc),, • S(Wc) 



(see (2.12 1 for notation). From these two equations one can deduce that satisfies 
Wc(x,'/)-X 9wc(x,??) dwc{x,v) 



1-77 



dx 



drj 



Then, using the connection between il(x,77) and wdXiV) (7.11), we see that the 
desired equation holds. □ 

There are certain boundary conditions satisfied by the complex slope J7(x,7?) 
which can be deduced from the study of local asymptotics (Theorem [2]) . Namely, 
we know that the limit shape has frozen facets outside the curve dV which is the 
same curve obtained as a frozen boundary in [KO' . Together with these boundary 
conditions, the complex Burgers equation ensures that the complex slope ri(Xj v) of 
the limiting ergodic translation invariant Gibbs measure coincides with that of the 
tangent plane to the limit shape of |CKP| . [KOj at (x, v)- That is, the 'hypothetical' 
limit shape for our model inferred from the asymptotics of correlation functions and, 
in particular, from the limiting density function Bn(0,0) = arg(f7)/7r (see Fig. [5]), 
coincides with the actual limit shape. 

7.7. Frozen boundary. Let us now discuss the frozen boundary curve dV. We 
obtain its explicit rational parametrization, and identify tangent points on it, thus 
proving Propositions |2.6| and |2.7| 



Proof of Proposition 2.6 When a point {XtV) G 2? C approaches the frozen 
boundary curve dT), the parameter Wc^XtV) — the critical point of the action 
S{w;x,v) (7.3) in the upper half plane — becomes real and merges with Wc(Xi'7)- 



This means that the frozen boundary can be characterized as a set (x, v) such that 
S'(u);x, 7?) has a double critical point. In fact, it is more convenient to take the 
double critical point Wj, itself as the parameter on dT>, and express x and 77 through 



it. We will now estabhsh the desired parametrization (2.11 )-(2.12 ). 
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Figure 13. Example of the frozen boundary curve with k — 12 
(mscribed in a 36-gon). 



The equations -^S{w;XtV) = :g^S{w,Xi''l) = for the double critical points 
of S{w;x,il) can be rewritten as (see (2.12) for notation): 



X 



n(wc 



1 



X 



X + 1 - 



S(w,). 



The first equation is just (7.8), and the second one is obtained by differentiating 
(7.12 ) with respect to w. The above two equations can be reduced to linear ones in 
X and rj, and solving them one arrives at (2.11 ). All other assertions of Proposition 
|2.6|are readily checked. 



□ 

From now on we will think of Wc, — oo < Wc < oo, as of the parameter of dV. In 
particular, the parametrization of dV by the double critical point ( 2.11 )-( 2.12) 
can be used to draw frozen boundary curves, see Fig. [4] and [T3| 

Let us turn to identification of tangent points on the 
The slope of the tangent vector {x{wc),ri{wc)) to dV is 



Proof of Proposition [£ 
frozen boundary curve, 
given by 



X(Wc 



- x(wc) 



77(Wc) 1 - T]{Wc) 



= ■ T(Wc). 



Indeed, this can be obtained by differentiating (2.11 ) with respect to Wf, and noting 

X(Wc) _ 



that n(wc) = n(wc)S(wc). One can then check that 

easy to verify using (2.11 )-(^2.12 ) that 

The rest of Proposition \~~ 
the tangent vector is related to w, 



n(wc) 



It is also 



7)(Wc) 1 — II(Wc) 

^^^^^ simplifies to the same expression. 



is not hard to check knowing how the direction of 
One should also use the fact that = . □ 

We conclude the discussion of the frozen boundary curve in this subsection by 
giving an illustration of how various parameters change along the curve. On Fig. [M] 
there is an example for k = 4. The horizontal axis is the parameter Wc. The straight 
line with slope 1 represents the double root of the equation £^S{w;x7'n) = 0. 
Two other roots are also shown. There are two points where Wc is a triple root, 
they correspond to turning points of the frozen boundary (one of them is marked as 
't'). The segments shown on the vertical line are [a^, bi]. The shaded region shows 
the segment [x(wc) + vi^c) — 1, x(wc)]- When Wc = ai,bi, or hi (see the marks 'a', 
'b\ and '/i'), the frozen boundary is tangent to a side of the polygon according to 
Proposition |2.7| 
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Figure 14. Various parameters along the frozen boundary curve. 
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Figure 15. Functions T{^Nc), ^{wc), and S'(^^(wc) along 91?. 



Fig. 15 shows the functions T{wc), ^(^c), and the third derivative S'^'^-'(Wc) :— 
^S{z; x(wc), r7(wc))|z=Wc as functions of the parameter Wc along the frozen bound- 
ary curve. The function T{\Nc) describes the direction of the tangent vector to the 
curve, and r2(wc) shows the type of the neighboring frozen facet outside the liquid 
region, see Proposition 2.7 and the discussion in ^2.5 (in particular. Fig. [?]). The 
third derivative changes its sign at every tangent and turning point where it is 
infinite or zero, respectively. 
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8. ASYMPTOTICS AT THE EDGE 



This section is devoted to studying the asymptotics of random lozenge tihngs on 
the edge of the hmit shape. See §2.5|for a prehminary discussion and illustrations. 



8.1. Limit regime. Fix a point {XiV) — {x{^c),v{'^c)) G dV which is neither a 
tangent nor a turning point. Here Wc is the parameter of the frozen boundary curve 
(Proposition |2 . 6[ ) . Let the points (xi, ni), . . . , {xs, Tig) be in the vicinity of the point 



N 



(i = 1, . . . ,s). 



In contrast to the bulk limit regime (Theorem [2|, at the edge to obtain nontrivial 
asymptotics, the differences Xi — Xj and rii — rij grow as TV — >■ oo. More precisely, 
we assume the following scaling: 

where 

x[ = 2-i/3(5(3))i/3(i _ _ n)-^(a, -Tf) ^ ' ' 

(recall that ^^^^ = ^(^^(wc) = S'(z; x(wc), r/(wc))U=wJ- Here n and cr^, i = 
l,...,s, are the new coordinates around {x,v)- Because of the factor T(wc) in 
(8.1) and due to (2.13), one sees that these new coordinates describe fluctuations of 
order N^^^ in tangent direction to dV and of order N^^^ in normal direction. The 
exact form of the scaling (8.2) is chosen to ensure that the limiting process would 
be the standard Airy process with kernel (2.15). 

8.2. Formulation of the result. Let K be the correlation kernel of our model 
(2.6). For r2(wc) > 0, let the scaled correlation kernel be defined by 

IC{xi,n,;x2,n2) -.^ N'/''\S'^''y2\'^''{l - 7j)\n\{l - n)-^K{xunuX2,n2), (8.3) 

and for il{\Nc) < wc set 

/C(a;i,ni;x2,n2) -.^ N'^'\S^''^ /2\'^'il - (8.4) 



l.Ti —cr 



'-xi=X2 lni=n2 



K{xi,ni;x2,n2)]. 



According to Proposition 2.7 this means that for the SV part of the frozen boundary 
(Fig. [t]) we are replacing particles (on Fig. |3| by holes and vice versa. For a general 
discussion of the particle- hole involution, e.g., see jBOOi Appendix 3]. In our model 
this involution is done because the neighboring frozen facet at the SV part of dV 
is densely packed with particles, and thus we must look at fluctuations of holes. 



Theorem 8.1. In the scaling (8.1)-(8.2), we have the convergence 
<lei[JC{xi,n.,;xj,nj)\l .j^^ det[^(Ti, (7^; r,-, ct-,)] ■ .,-^1, 



where A is the extended Airy kernel (2.15). 

The interpretation of this result in terms of fluctuations of nonintersecting paths 
is discussed in |2.5[ In the rest of this section prove Theorem |8.1| 
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8.3. Kernel K. Arguing as in the proof of Proposition 
regime 



7.2 



we see that in the 



K{xi,ni;x2,n2) 



l)-(8.2) the kernel K{xi,ni; X2,n2) (2.6| has the following asymptotics: 

P^%^ + ^x (8.5) 

1 exp(iV[S(w;;t,t)-^(^;t,t''^ 



dz 



dw- 



£(X-) c(oo) 



- z y/{w- x){w - x + l-r7)V(^- X){z - X + 1 - ??) 



The contours €{x~) ^-^id c{oo) are the same as described in Proposition 7.2 
It is convenient to conjugate the kernel K as 



K 



Aria/ ^2 "2% r.^ ^1 "1 ^^ 



,ni;x2,n2) 



(8.6) 



In the definition of the action S{\Nc', , ^^) ( 7.3 1 we must take cuts in logarithms 
in, say, the lower half plane because Wc S K. Of course, the new kernel K can also 



serve as a correlation kernel of our measure 



8.4. Expansion of the action at the edge. In the integral in K (8.6), there is 
an exponent of the difference of two expressions like 

7V(5(«;;f,f)-5(w,;f,f)). 

Let us expand such an expression in powers of N: 

Lemma 8.2. Let w = + w'N-^/^. For N{S{w; f , f ) - S{wa; f , f )) to be 
asymptotically bounded, it suffices to take 

x^xN + n'T{wc)N^^^ + x'N^'^, n = tjN + n'N^^^ 

for some x' ,n' G M. Then 



N{S{w- ^, ^) - 5(w,; ^, ^)) = i5(3)^.3 ^ 1 (1 _ ^)-2(^ _ nfn-^n' 



w ■ 



- (1 - r])-\i - n)^n'^x'w' + o(i). 

Proof. Observe that for partial derivatives at (yjcXiV) one has 



Denoting dx 
we obtain 



wxx — ^wwx^ ^wrjTj — ^wxv — ^wwrj- 

and taking the Taylor expansion at {\Nc'^x^v)^ 



n—Nrj 
N 



n(^S{wc + w'N 1/3; x + dx,V + dij) - S{sn^- x + dx,V + dr])^ 



(Swxdx + Swrfdr]) + ^Syj^^dx + Swrirjdxdri + ^Swijridr] 



We see that if dx, drj ^ const- N'^^"^ , and that also Swxdx + Sw-qdrj ~ const •iV~^/'^, 
then the above expression is asymptotically bounded. Then it is readily checked 
that the desired expansion holds. □ 



Lemma 8.2 justifies the scaling (8.1) of the coordinates {x,n), and also suggests 
scaling of variables in the double contour integral in the kernel K . 



ASYMPTOTICS OF RANDOM LOZENGE TILINGS VIA GELFAND-TSETLIN SCHEMES 41 



V 


u/ 


V 





Figure 16. Contours for the extended Airy kernel. 



8.5. Extended Airy kernel. We will use another formula for the extended Airy 

(8.7) 



kernel ^^(ti, cti; T2, 172) (2.15) which is written out in |BK1 §4.6]: 
-4(ri,cri;r2,cr2) 



1 



ri<T2 



^Att{t2 - Ti) 



1 



exp 



(gi - 0-2)^ 

4(t2 - Ti) 



ri)(cri +(72) + ^{t2 



(2^i) 



1 



1 



exp TlfTl - T2a-2 



((Tl - tI)u + (fT2 - 



TiU 



T2V 



3^ 



dudv 



Both contours come from infinity as straight lines, then meet at 0, and go to infinity 
in other directions. In u the contour goes f rom e~^^^^oo to to e^'^^^oo, and in v it 
is from e~^"^/^oo to to e^"^/^cx) (see Fig. 16). 

Comparing cubic polynomials in the exponent under the integral in (8.7) with 
the expansion of Lemma | 8.2[ we see that the parameters x'l 2 and n'l 2 should 
necessarily have the form (|8.2[ ) in order for our pre-limit kernel K to converge to 
the extended Airy kernel 



.7). Indeed, under 



1 - 



.2 ), we have: 



iV[S'(u'; 



N ■ 



^^)-^(z;t,t) + ^(Wc;t: 



)] 



'5'(wc, J ]\[ J ^y^j jv 

= -((Tl - )u + ((72 - t|)u - TlU^ + T2V'^ 

where u and v are the new variables, 

w ^\Nc + N-^S''^'^/2)-iu, z = Wc + N-^S''^'^/2)-h. (8.8) 
Thus, up to the factor e'^i'^i^'^^'^^-jT-j^+gTl (^.^^j^jgi^ s imp ly cor responds to another 
conjugation of the correlation kernel), we get in K (8.5)-(8.6) the same exponent 
as in the extended Airy kernel ( |8.7[ ). 
To make the new variables u and v 



.8) local around the double critical point 



Wc (as on Fig. 16), we need to transform the z and w contours so that they will 
cross at Wc and locally look as on Fig. [iTl A fter that, the leading contribution 
to the double contour integral for K ( |8.5[ )-( |8.6[ ) will come from a neighborhood of 
Wc, and we will see the desired convergence to the correlation functions of the Airy 
process (Theorem 8.1). 
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.6. Transforming the contours and comple ting the proof. Here we explain 



how the contours of integration for the kernel K (8.5)-(8.6) should be transformed. 
In order for the integration outside a neighborhood of w^, to be negligible, we require 
that on the contours ^S{w) < 5RS'(wc) < ^S{z) for z^w ^ Wc- In the process of 
moving the contours, certain residues may be picked, and as a result we obtain the 



desired convergence as in Theorem 8.1 



The rules for dragging the z and w contours are explained in the beginning of 
the proof of Proposition |7.8[ We consider three cases corresponding to various 
directions of the frozen boundary as indicated on Fig. [7] and in Proposition |2.7[ It 
is also informative to consult Fig. [M] and [15] 

8.6.1. Case VH: < il(wc) < 1 and bi < Wc < hi+i- In this case the double 
critical point Wc lies to the right of x because T{wc) > 0. Looking at curves 
{w: QS{w;XjV) — '^Slwcix^v)} ^-nd counting arguments (similarly to the proof 
of Proposition |7.8[ ), one can see that it is possible to transform the contours in a 
desired way as shown on Fig. [TSj 

For some configurations, it may be needed to drag the w contour through the z 



contour (as on Fig. 18 bottom). As a result, a residue at u; = z will be picked, but 
that residue will be integrated over a z contour not containing its poles which are 
inside [x 



V 



5 1 will arise. 



l,x]- Thus, no summand in addition to the one in 
For S'^'^) < 0, locally the contours around Wc will have directions that produce a 
minus sig n (se e Fig. 18 bottom). However, this sign is absorbed by the change of 
variables (8.8 1: recall that under the integral we have ^^2^. Due to the discussion 



8.5 



we have now established the desired convergence of the integral in K (8.5 1 



(8.6) to the integral in the extended Airy kernel (8.7) (with correction as in (8.9) 
below). 

Let us turn to the additional summand. In scaling 



1 )-( 8.2 ) , the fact that both 



conditions n2 < n-i and X2 < Xi are satisfied translates to the condition ti < T2. 
Then it can be shown using the Stirling app r oxim ation for the Gamma function 



(|7.6|) that the additional summand in K 



51 

Ma-v) N(s( 



6|) behaves as 



X2 "2, 

' N ' N ' 



)). 
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Figure 18. Case VH (j |8.6.l| ). S'^^) > O (top) and S^^'> < (bot- 
tom). Transformed contours (left) and curves where '^S{w; x, ^7) = 
5S'(wc;x,»7) (right). 



(Xi - X2 + l)ni--;t2- 

(ni - 71.2 - 1)! 



CXP (-5(^1 - Cr2 - tI + t|) V(t2 - Ti 

^/47r(r2 - ri) 



Summarizing, we conclude that in the case VH the rescaled kernel (8.3) converges 
to the extended Airy kernel as follows: 



/C(xi,ni;x2,n2) ^ e-^^^^+^^'^'+i^'-i^^ A{ti, ai;T2, a2). 



(8.9) 



Because the exponent factor is just a conjugation not affecting the determinant, 
this implies Theorem |8.1| in this case. 



8.6.2. Case HS: il{wc) > 1 and hi < Wc < 



(8.1)- 



First, observe that the scaling 
.2 ) implies that the two conditions n2 < ni and X2 < xi cannot hold 

Thus, there is no additional summand in K 



simultaneously because T[wc) 
coming from 



< 0. 



.5). 



Since in this case T(wc) < —1, the double critical point lies to the left of 
X + 77 — 1. We transform the contours as on Fig. 19 (one can also draw curves with 
95(w;;x,»7) = '^S{wc;XtV) as on Fig. 18 but we omit this): 

• For S^^^ < 0, we replace the counter-clockwise z contour coming from €{x2) by a 
clockwise contour corresponding to €'{x2) (see Lemmas 6.2 and 6.3 for notation). 
This leads to appearance of the residue ( |6.11 1 . The change of sign in the integral 
due to local directions of the contours around will be again compensated by 
the change of variables (8.8). 
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— ^ 1 







Figure 19. Case HS (^8.6.2). Transformed contours for S^^^ < 
(left) and S^^^ > (right). 



For S^^^ > 0, the w contour is dragged inside the z one, and the resulting residue 
at w = z gives ( 6.10 ) after the integration. Alternatively, one could also transform 



the contours in the same way as in the S'*^'^-' < case. 
We see that either way, a transfor mat ion of con to urs results in the additional sum- 

B.2), the condition ni > n2 turns 



mand l„i>,: 



(2;i-a;2 + l)„ 



5). By 



into Ti < T2. It can be shown that for our case fi(Wc) > 1, the rescaled additional 
summand 



(ni - ^2 - 1)! 

converges to -l^i<^2exp (-Kdi - cr2 - + T|)^/(r2 - n)) /\/47r(r2 - n). This 

fact and the transformed contours ensure the desired convergence (8.9), which im- 
plies that Theorem |8. 1| holds in the HS case as well. 



.6.3. 



Case SV: ri(wc) < a nd a j < vj^ < hi. As in j 8.6.21 th er e is no additional 



1H8.2) 



summand in K coming from (8.5) because T{\Nc) < 0, see 

Because — 1 < T{^c) < and Ui < Wc < h, we see that the double critical point 
Wc belongs to the intersection of the segments (x + '7 — 1, x) ^1 (^^ij ^O- We transform 
the contours as shown on Fig. [20] (we also omit the curves where one can also draw 
curves with '^S{'w; x, 77) — 35(wc; Xi v))- As before, for 5*^3) < q, the change of sign 
due to directions of the contours is compensated by the change of varia bl es |8.8 [). 
There would also be another minus sign in front of the integral i n K (8.5 )-(8.6 ) due 
to the factor {{w - 



■x)iw-X+^~v)iz-x){z-X + '^-v)y 



5). This justifies 



the need of the particle-hole involution for il(wc) < (se e (|8.4[ )) at the level of 
formulas. So, scaling the kernel K after the involution as in ( |8.4[ ), we conclude that 
for the double contour integral part of the kernel, the desired convergence holds. 

While transforming the contours, we drag the w contour through or inside the 
z one. This results in the residue at w = z. It will be integrated as in Lemma 



6.2 



and this will give an additional summand 1 



(rii-n2)x 



particle-hole involution 



X2<xi (2;i-a;2)! 



laji — X2 Ini - 



1 



X2<Xi ' 



in 



4), we see that 

— "•2)2:1-3:2 



= -1 



[Xi - X2) 



Xi >X2 ' 



(ni - n2)xi-x2 
(xi - X2)l 



Applying the 



^xi—X2^ni^n2 • 



For ri(wc) < 0, the scaling (8.1 1-(8.2 ) implies that the conditions Xi — X2 and rii ^ 



n2 cannot hold simultaneously. One can then check that the following convergence 
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Figure 20. Case SV (j |8.6.3[ ). Transformed contours for S'^^) > 
(top; there are variations depending on relative positions of [a^, 5^] 
and [x + ^ Ij x]) ^-nd 5*'^^ < (bottom). 



holds: 
-1 



. Arl/3| C(3) /o|l/3 1^1(1 -^) ,iv(s(w.;f,f)-S(w.;f,^)) K-^2k-.. 

exp -02- tI + t|) V(t2 - Tl) 



^ -1 



Tl<T2 " 



\/47r(r2 - Tl) 

Thus, Theorem |8.1|is established in the last remaining SV case. 
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